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ABSTRACT 



This research considers a new realization of digital filters suitable for VLSI 
implementation. The method involves delta modulation which provides analog-to- 
digital (binary) conversion. The output of a linear system is the convolution of the 
input and the system impulse response. This new digital filter requires that both 
the input and the impulse response be first converted to bit streams using delta 
modulation. These bit streams axe then convolved. The result is an analog voltage 
which approximates the convolution of the analog functions. 

Direct convolution of the bit streams is difficult to realize with electrical cir- 
cuits. A greatly simplified system with equivalent performance is a result of this 
research. This is called the reduced delta convolution (RDC) system (digital filter). 
The performance of the RDC system when used as a convolver and as a correlator 
is analyzed and verified by computer simulation. Analyses of the effects of self 
noise and external noise are included. Conclusions axe that the RDC system has 
considerable potential as a digital filter when using integrated circuits. Realization 
requires considerably fewer components and simpler connections than other digital 
filters. A reason is that there are no multipliers required in the RDC system. The 
RDC system requires no synchronization, operates in real time and is easily pro- 
grammed. Further, the RDC system has noise performance which is better than 
predicted by ordinaxy filter theory. 
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I. INTRODUCTION 



This research is concerned with a new realization of digital filters. The method 
used involves delta modulation. 

We define and discuss digital filters in Chapter Two. 

In Chapter Three we briefly explain delta modulation which converts an analog 
signal to a digital (binary) signed. We then introduce a new method of convolution 
which is called delta convolution in this report. The two signals to be convolved 
are first converted to binary form using delta modulation. Time domain analysis 
is used to derive equations of interest and to define the required hardware. An- 
alytical examples are presented to show that the results obtained using the delta 
convolution method and conventional methods are similar. We conclude that direct 
implementation of the delta convolution method results in a complicated system 
(large number of gates and interconnections). 

A system which involves two feedback circuits to account for the previous 
values involved in delta modulation and convolution is developed in Chapter Four. 
We call this the reduced delta convolution (RDC) system. The result is a simple 
hardware realization. We define (block diagram) a system which functions as a 
convolver (linear filter) and a similar one which is a correlator. 

In Chapter Five we introduce a possible circuit to implement the RDC system. 
A comparison between the needed hardware and complexity of the RDC system 
and a conventional digital matched filter (DMF) is included. The result is that the 
RDC system has a simple layout for IC realization and requires less than 10 percent 
of the transistors needed for a typical DMF having comparable performance. 
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Simulation results when the RDC system is used as a convolver (lowpass filter) 
and a correlator (matched filter) are presented in Chapter Six. The simulation 
results include the transfer function of a RDC lowpass filter, the response of the 
RDC lowpass filter to two signals, one in the pass band and the other outside 
the pass band, and the output of a RDC lowpass filter with a square wave input. 
Included also axe responses of matched filters to a raised cosine pulse and to a chirp 
signal. 

Chapter Seven documents the noise characteristics of the RDC system. Op- 
timum values of important RDC system parameters are derived. Analyses of self 
noise and input noise are presented. A conclusion made is that the RDC system 
has noise performance which is better than predicted by ordinary filter theory. 
Simulation results for the chirp matched filter and lowpass filter when used in the 
presence of additive input noise are shown. Again, the noise performance is better 
than predicted by theory. A reason for this is detailed. This chapter also includes a 
modification to the RDC method called variable step size RDC. Simulation results 
for the variable step size RDC lowpass filter and notch filter are presented. 

The results of this research indicate that the VLSI compatible RDC system 
may be an attractive means of realizing digital filters. The system requires no 
synchronization, has no multipliers and is easily programmable. Other features 
include real time operation (no off-line processing) and no limitations on the time 
duration of the signal to be represented in matched filter applications. 

When compared with other digital filter realizations, the RDC system appears 
to require considerably fewer components and simpler connections. Finally, because 
of the properties of delta modulation, the RDC system achieves noise performance 
which surpasses that predicted by ordinary filter theory. 
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II. BACKGROUND 



This research is concerned with the realization of linear systems (filters) using 
digital technology. Common realizations of digital filters require a large number 
of transistors and interconnections. This report describes a new type of digital 
filter and documents its advantages and performance. We begin this description 
by introducing linear systems. 

A. THEORY OF LINEAR FILTERS 

A continuous-time or analog signal may be described by a function of time, 
say x{t). The response of a linear time-invariant system to x(t ) is given by 

OO OO 

y(t) = / h(r) ■ x(t — t) ■ dr = J x(A) • h(t — A) • d\ (2-1) 

— OO — OO 

where h{t) is the response of the system at time t to a unit impulse applied at time 
t = 0. See Fig. 2.1. The function h(t) is, then, used to represent the linear system. 
Linear systems are used as convolvers to realize filters (lowpass, bandpass, notch) 
and as correlators (matched filters (MF)). 

The right side of Equation 2.1 is commonly called the convolution of two 
functions x(t) and h(t ) and denoted as x(t) * h(t). From transform theory 

Y(f) = X(f) ■ H(f) (2.2) 

where and H(f) are respectively the Fourier transforms of y(t),x(t ), 

and h(t). 
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A discrete-time signal is a sequence of numbers {x(? • T s )}, where the integer 
i may vary over a finite or an infinite range, and where T s is the sampling interval. 

A linear time-invariant discrete-time system can be described by the input- 
output relationship 

OO 

y(n ■ T s ) = T S ■ x ( m -Ts) • h((n - m) ■ T s ) 

m = — oo 

~ - OO — 

= ?*■ *((n-k)-Z)-h(k-T a ) (2.3) 

— oo 

where x(n ■ T s ) and y(n • T s ) are the input and output signals, and h(n ■ T s ) is the 
impulse response of the system. That is, h(n ■ T 3 ) is the response of the system at 
n • T s due to a unit sample (unit impulse sample) applied at t — 0. See Fig. 2.2. 

The right side of Equation 2.3 is the convolution sum of the two sequences 
{x(n-Tj)} and {h(n-T 3 )}. When the sequence {h(n-T 3 )} has only a finite number 
of non- zero terms, we say the system has a finite impulse response (FIR). Otherwise, 
the system is said to possess an infinite impulse response (HR). If h(n -T s ) = 0 for 
n < 0, we say that the system is causal or physically realizable. 

Using ^-transforms 

Y(z) = H(z)-X( z ) (2.4) 

where Y(z ), H(z ), and X(z) are respectively the z- transforms of y(n ■ T 3 ), h(n ■ T 3 ) 
and x{n • T a ). 

A large class of linear time-invariant discrete-time systems can also be de- 
scribed by the linear constant coefficient difference equation 

N L 

y(n-T s ) = y>- -y{{n- k) ■ T 3 ) + ^ b k ■ x((n-k)-T 3 ) (2.5) 

k=l k = 0 

where it is possible to .convert Equation (2.5) to an equation of the form of Equation 
(2.3) (Kirk and Strum [1]). 
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output 
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Fig. 2.1. Representation of linear time-invariant continuous- time system. 



input 
x(n .T s ) 






h( n .T s ) 



output 

► 



y( n -T s ) 



Fig. 2.2. Representation of linear time-invariant discrete-time system. 
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In Equation (2.5), if a*, = 0 for all k, then 

L 

y(n -T a )= k) ■ T s ) (2.6) 

k= o 

So, the output of the system depends only upon the current and previous values 
of the input. This system is called non-recursive. 

From these results, we can now consider ways of realizing discrete-time sys- 
tems. These realizations are called digital filters. This research is concerned with 
a new realization of a digital filter. 



B. VARIOUS REALIZATIONS OF LINEAR FILTERS 



The unit sample response for the non-recursive system described in Equation 

(2.6), is given by 

L 

h(n • T 3 ) = Y / bkS(n-k). (2.7) 

k = 0 

This system is known as a finite-impulse response (FIR) system because we have 
a finite number of terms in the unit sample response (Bird [2]). 

From Equation (2.7), a FIR digital filter has a transfer function 



L 

H{z) = Y2, bkz ~ k 

k= 0 

which is often implemented as a tapped-delay line and is sometimes called a 
transversal filter. (See Fig. 2.3 and Rabiner [3]). 

It is important to realize that, while a non-recursive filter must have a FIR, 
a FIR filter need not be implemented non-recursively. A recursive implementation 
of a FIR filter can be obtained by introducing poles in the transfer function and 
then canceling them with extra zeroes (Bird [2]). A new recursive implementation 
of a FIR filter, which is called the reduced delta convolution (RDC) method, is 
introduced in Chapter Four. 
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Fig. 2.3. Realization of transversal filter. 



A transversal filter can have an arbitrary impulse response of finite time dura- 
tion and therefore can be used to implement any linear filter. (Any system whose 
output is linearly related to the input is a linear filter.) In this sense, a transversal 
filter can be thought of as a fundamental building block of linear systems. 

In Equation (2.5), if a* ^ 0 for one value at least, then we have an HR filter. 

In classical signal processing most of the filters are HR systems. These filters 
correspond to direct implementation of the filter’s rational polynomial transfer 
function. 

In modern signal processing most filters are FIR systems. These filters corre- 
spond to a direct implementation of the filter impulse response approximated with 
a finite number of terms. These filters often arise as matched filters for radar and 
sonar signals (Kailath [4] and Rabiner [3]). 

FIR linear filters have several important properties which make them attrac- 
tive for digital signal-processing applications. Among these features are simple 
design, linear phase and the absence of any stability problems which may occur 
in HR filters. On the other hand, for the same filter, non-recursive technique uses 
more taps and multipliers (Hamming [5], Young [6], McClellan [7], and Whalen 
[8])- 

' Transversal filters are the popular choice of digital filter when using VLSI. We 
consider next various circuit realizations of transversal filters. 

C. TRANSVERSAL FILTER REALIZATION 

1. Block Diagram 

a. Transversal Filter as a Linear Filter 

To implement a linear filter using transversal filters, it is sufficient to 
choose the tap weights appropriately. From Equations (2.6) and (2.7), if the desired 
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impulse response is known, then the transversal filter tap weights are simply the 
sample values of the filter impulse response. See Fig. 2.4. 

b. Transversal Filter as a Matched Filter/Correlator 

Pulse detection in noise differs from analog signal transmission in two 
major respects. First, we are usually concerned with determining the presence 
or absence of a pulse. Second, we often know the pulse shape in advance, but 
not its amplitude or arrival time. The pulse-shape information makes it possible to 
design optimum receiving filters for detecting pulses buried in noise having a known 
spectral density function. Such optimum filters (correlators) are called matched 
filters. 



Let x(t),0 < t < T d , be a signal to which a filter is to be matched. 
The impulse response of the matched filter is defined as (Turin [9] and Turin [10]) 



( b • x(T d - t); 0 < t < T d 

h(t) = < 

1 0; elsewhere 

where the gain factor b is arbitrary and is henceforth taken as unity and T d is the 
pulse time duration. The name matched filter comes from the fact that h(t) has 
the same shape as the pulse x(t ) with time reversed and shifted by T d seconds. 
Fig. 2.5 is a diagram showing the application of a matched filter. 

The ratio of peak signal power S p to average noise power N 0 at the 
output of the MF is given by 



k mT ‘- B r t 



where Si and N{ are average input signal power and noise power and B is signal 
bandwidth. Then 

SNRI = = T d B (2.8) 

where SNRI is the signal-to-noise ratio improvement factor (Horrigan [11]). 
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Fig. 2.4. Diagram of a transversal filter as a linear filter. 
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matched filter 




Fig. 2.5. Diagram indicating the application of a matched filter correlator. 
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This important result shows that the filter SNRI is determined only by 
the time duration of the pulse and its bandwidth and not by the pulse ‘shape’. 

Fig. 2.6 shows a transversal filter implementation of a MF, where we can 
see that the values of the weights are exactly the same as that given for linear 
filters (Fig. 2.5) except that the positions of the weights are reversed. 

2. Circuit Realization 

In Figures 2.4 and 2.6 the delay needed for the analog signal can be im- 
plemented using the following conventional techniques: 

• Charge- Transfer devices (CTD’s), which include both charge-coupled devices 
(CCD’s) [12 - 14] and bucket-brigade devices (BBD’s) [15, 16], can be used to 
delay the analog signals (Butler [17]). When CTD’s are used to delay analog 
signals, the signal to be delayed is first sampled at a rate greater than twice 
the largest significant frequency component of the signal. The analog samples 
are then clocked through the CTD shift register. CTD’s can achieve hundreds 
of milliseconds of delay (Buss [18]). 

• Acoustic delay lines are used as alternatives to CTD’s for analog time delay, 
but for small values of time delays (< 20^zsec (Buss [18])). 

• Digital shift registers (SR) preceded by analog/digital conversion (ADC) and 
followed by digital/analog conversion (DAC) can also be used to delay signals 
as shown in Fig. 2.7. 

a. Transversal Filters Using CTD’s 

In order to make a CTD transversal filter, it is necessary to nonde- 
structive^ sample the delay line and to perform the weighted summation. The 
sampling is achieved in different ways depending on whether CCD’s or BBD’s are 
used (Buss [18]). 
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Fig. 2.6. Diagram of a transversal filter as a matched filter. 




Fig. 2.7. Diagram showing realization of time delay of analog signals 

using digital circuits. 
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Generally, there are two approaches to implementing programmable 
transversal filters using CTD’s. 

The first is the analog/ analog system (Bosshart [19], and Denyer [20]). 
Fig. 2.8 shows a block diagram of an analog/analog CTD transversal filter where 
analog stores for both input signal and impulse response (reference store) are used 
in conjunction with analog multipliers. 

The main limitation is that since the weighting coefficients are stored 
in the form of charge, they decay due to thermal leakage and must be refreshed 
every 10 to 100 ms. The tap weight resolution is limited by the fixed pattern 
noise associated with the MOS analog multipliers at each tap. In addition, the tap 
weight data that are stored in off-chip digital memory must be converted to an 
analog signal before it can be transferred to the on-chip analog reference store. 

The other technique is the digital/analog approach (Tiemann [21] and 
Tanak [22]) where the weighting coefficients are represented in digital form. See 
Fig. 2.9. This requires a multiplying digital-to-analog. converter (MDAC) at every 
filter point. The MDAC uses MOS transistors that route the output of floating 
gate taps to either positive or negative summing buses. As a result of this, a fixed 
pattern noise results due to variations in the tapping transistor characteristics along 
the length of the array (Gandolfo [23]). 

We can summarize the main limitations of transversal filters realized 
using CTD’s as follows: 

• The time duration of signals that can be processed using CTD matched 
filters is ultimately limited by the storage time of the devices (i.e., the time it 
takes a stored charge to be lost due to leakage). Typically the storage time is 
of the order of 1 second (Buss [18]). 
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OUTPUT 



D is a delay stage 



(a) 




Fig. 2.8. Diagram of the analog/analog method of implementing CTD 
transversal filter, (a) System, (b) System impulse response. 
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prefilter sampler 




Fig. 2.9. Diagram of the digital/analog method of implementing a CTD 

transversal filter. 
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• The filter length (i.e., the number of delay stages L) is ultimately limited by 
charge transfer inefficiency. Calculations indicate that 



L • e < 2 



(2.9) 



where L is the number of delay stages and e is the charge transfer inefficiency 
per point where 



and 



e = 10 3 for early versions of CCD/MOS 



e = 0.4 x (10) 3 for recent versions 



(Creasey [24]). 

• The Nyquist sampling theorem requires that a signal having a bandwidth B 
be sampled at a frequency greater than 2 B. Combining this requirement with 
Equation (2.9) gives the following limitation on the Td • B product (which is a 
measure of pulse compression ratio or processing gain) of signals that can be 
processed using CTD filters. 

Td ■ B < 1/e. (2.10) 

The limitation imposed by Equation (2.10) can be compensated for by selecting 
the weighting coefficients to invert the dispersion due to the imperfect charge 
transfer. But the dependence of charge transfer efficiency on the signal ampli- 
tude makes it impossible to exactly invert this dispersion at all signal levels. 
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• The signal bandwidth is limited to less than half the maximum clock frequency 
of the filter. 

o For CCD’s this limitation is « 20 MHz. 
o For BBD’s this limitation is « a few Megahertz. 

• Another limitation on the CTD transversal filter is the accuracy with which 
the weighting coefficients can be determined. For example the transversal 
filters using MOSFET analog multipliers suffer from: 

o Poor accuracy 
o Lack of long-term stability 

o Drifts that axe found in most MOSFET multipliers 
The weighting coefficients can also be realized using a split-electrode structure. 
However, because the weighting coefficients are fixed by the split-electrode 
structure, such devices are only suitable for applications where a fixed filter 
response is required (Haken [25]). 

b. Transversal Filters Using Digital Delay 

It is possible to use clocked shift registers to provide delay. Then in 
Fig. 2.10 the sampled input x(n • T s ) is quantized into a number of bits, typically 
12 to 16 bits for audio and 6 to 8 bits for video signals. 

The digital signals (code words) recirculate in virtual delay lines com- 
posed of random access memory (RAM) and are multiplied by the weights which are 
stored in programmable read only memory (PROM). Every sampling time period 
a new code word enters the RAM where the oldest one is discarded. Using a single 
multiplier, the multiplication rate has to be L • f 9 , where L is the filter length and 
f 9 is the sampling frequency. The large number of multiplications required each 
second is the main disadvantage of this scheme (Kailath [4]). 
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Fig. 2.10. Diagram of a digital transversal filter which uses a single multiplier. 



Fig. 2.11 shows an equivalent scheme which uses L multipliers and 
where m.\ and m 2 are the code word lengths of the PCM encoder and micropro- 
cessor system respectively. In this scheme the multiplication rate is reduced to 
f s . The main disadvantage of this scheme is the increased cost compared with an 
equivalent CTD implementation (Terrell [26]). 

The next chapter introduces a new approach to the realization of transversal 
filters. The circuit that results requires no voltage multipliers which results in 
a considerably reduced number of components (transistors) when compared with 
CTD and digital delay versions. 
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III. DELTA CONVOLUTION 



A. INTRODUCTION 

This chapter introduces a new method of convolution which is called DELTA 
CONVOLUTION in this report. The two signals to be convolved are delta modulated 

B. DELTA MODULATION 

Delta modulation (DM) represents an analog waveform with a binary sequence 
(bit stream). This analog-to-digital conversion (ADC) is quite different from pulse 
code modulation (PCM). PCM involves code words and requires somewhat com- 
plex encoders and decoders. Recovering the analog waveform from the code words 
requires synchronization. By contrast, DM involves no code words, requires no 
timing signals, and uses simple hardware (Roden [27]). 

It has been found that analog signals such as speech and video signals gen- 
erally have a considerable amount of redundancy; that is, there is a significant 
correlation between successive samples when these signals are sampled at a rate 
higher than the Nyquist rate. The redundancy in these analog signals makes it' 
possible to predict a sample value from preceding sample values and to trans- 
mit the difference between the actual sample value and the predicted sample 
value estimated from the past samples. This results in a technique called DIF- 
FERENCE ENCODING. DM is one of the simplest forms of difference encoding 
(Ha (281). 
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The idealized DM codec (coder/decoder) is shown in Fig. 3.1. The band 
limited output of the LPF is compared by subtraction with the stepwise approx- 
imation x s (t). The difference x A is passed through a hard limiter whose output 
is 



(+ l \ : f ( x b(t) > x s(t)\ 

v-i y 11 wo <*.(*)/• 

The resulting binary output x(t) is applied through an ideal integrator (or accu- 
mulator) with a feedback gain factor G to produce x s (t) where 

{ Ax when using an accumulator 

Y~ when using an integrator; T s is the sampling interval ' ' ' 



and where 



x 3 



(2n + l 

V 2 




Ax • E * (^p-Ts) ; x (2^7,) e {— i, + 1 } 

i= 0 

when using an accumulator (3.2) 
Y~ • J n x(t) • dt when using an integrator 

^ — OO 



where 



% 3 



f 2n -(- 1 
\~~2~ 




; n ■ T s < t < (n + 1) • T s 



Vn 



The decoder for a DM waveform is simply a staircase generator (accumulator) 
or an integrator. If a “1” is received the staircase increments positively. If a “0” 
is received the staircase increments negatively. The accumulator (or integrator) is 
followed by an amplifier having an amplification factor G. The amplifier may be 
followed by a low pass filter (LPF) to smooth the staircase output into a continuous 
function. A set of waveforms associated with the DM are shown in Fig. 3.2 which 
also defines Ax of Equation (3.1). 

DM systems are subject to two types of quantization distortion which are gen- 
erally referred to as granular quantization noise and slope overload noise. Since 
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Fig. 3.1. Block diagram of the delta modulation codec, 
(a) Encoder, (b) Decoder. 
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Fig. 3.2. Delta modulation waveforms. 
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x 3 (t ) cannot change by more than Ax units in T 3 seconds, ^ is the highest input 
signal rate-of-change that the DM codec can follow. We call ^r- the DM slope- 
following capacity and denote it by x' Q . When \^\ exceeds this quantity, slope 
overload occurs and gives rise to the kind of error shown in Fig. 3.3. The granu- 
lar noise arises because the DM signal is a discrete- amplitude representation of a 
continuous amplitude process (Goodman [29] and Greenstein [30]). Consequently, 
granular noise is always present. Only granular noise occurs when the input is 
changing very slowly (~ constant). 

The key to effective use of delta modulation is the intelligent choice of two 
parameters: (1) the step size Ax, and (2) the sampling interval T s . Obviously, 
granular noise is reduced by decreasing Ax, but at the expense of a reduced slope- 
following capacity, and, hence, greater slope overload noise. In Chapter 7 we 
introduce a procedure to choose the optimum step size and sampling interval such 
that the total noise obtained from the granular noise and the slope overload noise 
is minimized. 

C. DELTA CONVOLUTION 

The subject of this research is the representation of any linear filter using 
DM. Time domain analysis is used to derive equations of interest and to define 
the required hardware. In this chapter we pursue the DM equivalent of a linear 
filter by direct operation on the convolution integral. This results in a cumbersome 
expression and complicated hardware. However, working from this first result, it 
is possible to apply two successive modifications which result in a quite simple 
equation and a simple hardware structure. These modifications are presented in 
Chapter 4. 
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Fig. 3.3. Waveforms showing slope overload and granular noise. 
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Assume that h(t) and x(<) are two analog signals where 



h(t) = x(t) = 0; 



t < 0 



Let h(t) and x(t) be the delta modulated versions of h(t) and x(t), where h 3 (t ) 
and x s (t) are the staircase representations of h(t) and x(t) such that we have no 
slope overload (Fig. 3.2). Assume also that 



h{ = i— bit of h(t) = h(t) for i ■ T s < t < (i + 1)T S 



and 



and 



X{ = i— bit of x(t) = x(t) for i • T a <<<(? + 1) • T a 

We let h t and x,- G {— 1,+1} for all i. Therefore, 

j 

h 3 (t ) = Ah • ^2 hi', j • T 3 < t < (j + 1) • T 3 

i= 0 

j 

x a (t) = Ax - ^2 Xi ; j • T, < t < (j + 1) ■ T. 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



t=0 



When the sampling rate and step sizes are chosen properly, then, as seen in Fig. 3.2, 



x(t) « x 3 (t); j • T a <t<(j + 1) • T 3 . 



Therefore, from Equations (3.6) and (3.7) we will have 



x(t) « x 3 (t) = Ax • ^2 j • T a < t < (j + 1) • T s 



(3.7) 



(3.8) 



i=0 



where 



j = integer part of 

In the same manner, 

j 

h(t) R3 h a (t) = A h-'^hi] j -T 3 < t < (j + 1) • T a 



(3.9) 



i=0 
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We now have expressions for x(t) and h{t) which involve the bits of x(t ) and h(t). 

Now, a linear system (filter) is characterized by its impulse response h(t). 
Given an input voltage x(t) and impulse response h(t), the output voltage y(t ) can 
be written as the convolution of x(t ) and h(t ). That is, 



y(t) = h{t ) * x 




— oo 



h(r ) • x(t — r) ■ dr 



A conclusion is that in the time domain, linear filters are completely and uniquely 
identified by their impulse response h(t). 

The preceding development suggests that x(t) « x s (t ) and h(t ) ss h s {t ) and 
so we can write 



y{t) K, h s (t) * x s (t) = J h s (r) ■ x s (t - t) ■ dr 



(3.10) 



Let k = integer part of . Then from Equation (3.9) 



hsi^) — A/i • ^ ^ /ip, 

p = o 

Similarly, from Equation (3.8) 

m 

x s (t — r) = Ax • X{- 



k • T s < t < (k + 1) • T s 



where m = integer part of 



i=0 



t—T 

T, 



But j = integer part of 



T, 



(3.11) 



m-T 3 <t — T<(m + 1)T S (3.12) 



and so 



t = j ■ T s + E i| €i is positive and less than T s . 



But k= integer part of 



T 

T, 



and so 



r = k • T s + € 2 ; €2 is positive and less than T s . 



(3.13) 



(3.14) 
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Prom Equations (3.13) and (3.14) 



integer part of 



t — r 



= integer part of 



0 — k) • T s + (€i — € 2 ) 



T s 



because €1 and €2 are positive and less than T s . So, 



the integer part of 



t — T 



= j — k 



(3.15) 



From Equations (3.12) and (3.15) 



j-k 



,(t - r) = Ai ■ ^ x,-; (j - k) -T a < t — t < (j — k + 1 ) • T a (3.16) 



«=o 



From Equations (3.10), (3.11) and (3.16) 



y(t) 



f 



A h-J2 h t 

P = 0 
k 



*f A h-Y,h p 

J L p = 0 . 



j-k - 

Ax • ^ Xi 

t=0 

j-k 

Ax • x,- 

t'=0 



dr 



• dr 



(3.17) 



where L\ is the time duration of non-zero values of h(t). A duration of L x seconds 
corresponds to n x = bits of h(t). 

A discrete version of y{t ) is formed by replacing the integral in Equation (3.17) 
with another summation to obtain 

k 



"i-l ( 

»(0»#([>'+i]-r,) = Ai-AA.r.. Y, \ 

k=Q l lp=0 
n\ — 1 ( k 

y(t) ~ y(j ■ Ts) = Ax • A/i • T s ■ h f 



j-k 

Li=0 



fc=0 t Lp=0 



j-k - 1 

Z X *' 

L i=0 



j = 0 , 1 , 2 ,.. 



; = 1 , 2 , 3 , 



= 0 for j < 0 (3.18) 

We write y(t) « y(j • T s ) = C x • jd(r)Ax(; • T s - r)| ; j = 1,2,3, . . . 



(3.19) 



= 0 



i < 0 
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where 



C\ = Ax • A h ■ T s 



and 



ni-l ( .. 

h(r)Ax (j • T s - r) = ^ l ^ h i 



k= 0 V Lp=0 



[i-k-i • 

Y x p 

- p= o 



where h p and x p € {—1, +1} . 

We call h(r)Ax(j • T, — r) the delta convolution (normalized) at time j -T 3 . In 
this last form, y(t ) can be represented as the sum of appropriate products of ±1. 
This result provided the initial motivation for this research because the product 
(±1) • (±1) can be formed with XNOR gates. This is a considerable hardware 
advantage compared with usual digital multipliers (and usual digital filters) which 
have to form the product of PCM code words. 

As an example of delta convolution, we convolve h(t) with x(t) where h(t) and 
x(t) are given in Fig. 3.4. Applying the conventional methods of convolution the 
result y cvi (t) is given in Fig. 3.5 when we let T s = 1. 

Assume that Ax = Ah = 1. Then, the delta modulation versions of h(t) and 
x(t) are given in Fig. 3.6. Applying the delta. convolution Equation (3.18), 



y(o) = o. 



The delta convolution when j = 1 (Fig. 3.7) is given by Equation (3.18) as 

y(T 3 ) = C, • {h(r)Ax(T 3 - r)} 

= C i * {/io * £o } 

= C, .{11} 

— C i = 1 with our assumptions. 
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Fig. 3.4. Waveforms, (a) Impulse response, (b) Input signal. 
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Fig. 3.5. The result of convolution using the conventional methods. 
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Fig. 3.6. Delta modulation versions, (a) Impulse response, (b) Input signal. 
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Fig. 3.7. The needed products for delta convolution at time T s . 
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The delta convolution when j = 2 (Fig. 3.8) is given by Equation (3.18) as, 

y(2T a ) = C 1 -[h(r)M(2T a -r)} 

= C\ • {}iq(2xq + £ 1 ) + h\ lo } 

= C\ • {1(2 + 1) 4- 1(1)} 

= C\ • (4) = 4. 

The delta convolution when j = 3 (Fig. 3.9) is given by Equation (3.18) as 

y(3T s ) = C x ■ [h(r)Ax(3T 3 - r)} 

= C\ • {/io(3£o 4- 2ii + £ 2 ) + ^i(2£o + Xi) + h,2XQ } 

= C\ • {1(3 + 2 + 1) -j- 1(2 + 1) + 1(1)} 

= Cl • 10 = 10. 

In the same way we can find the result of delta convolution y{t) for all values 
of j • T s (Fig. 3.10). It is clear from Figures 3.5 and 3.10 that simulated results of 
analog convolution and delta convolution are essentially the same for this example. 
Note that the maximum value of the convolution is 344.0 for both simulations. 
Theoretically, the maximum value of the convolution is 341.33. 

It is possible to expand Equation (3.18) and then create a circuit which forms 
y(j-T 3 ) for any h(t) and x(t). This is done in Appendix A. The result is a circuit 
having an unacceptably large number of XNOR gates and a large memory. We 
present in Chapter 4 a much simplified form of Equation (3.18) which eliminates 
the considerable redundancy of mathematical operations inherent in that equation. 
The practical result is a simple hardware realization of a circuit to form y(j ■ 
Ts). 
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Fig. 3.8. The needed products for delta convolution at time 2 T s . 



55 




Fig. 3.9. The needed products for delta convolution at time 3 T a . 
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Fig. 3.10. The result of convolution using the delta convolution method. 



IV. REDUCED DELTA CONVOLUTION 



Implementing Equation (3.18) directly results in a large number of gates and 
interconnections as noted in Appendix A. The following reasoning suggests a sim- 
pler form may be possible. The delta modulation representation of a function ( x(t ) 
or h(t)) at any time t involves all previous values (history) of the function. The 
feedback portion of the delta modulator accounts for this history. Also, the result 
at any time t of convolving two functions depends on the history of both functions. 
The integral portion of the convolution equation accounts for this history. 

When using Equation (3.18) to determine the delta convolution result at time 
t = < 2 ,^ 3 , etc., we essentially involve (calculate) repeatedly the previous history 
of both x(t) and h(t ) at *2,^3, etc. This implies considerable redundancy of cal- 
culation (hardware) which increases with time. By carrying the history (previous 
calculations) forward, to the next calculation, the required hardware should be re- 
duced. Feedback can be used for this purpose. Indeed, use of one feedback circuit 
does reduce the complexity as shown in Appendix B. In this chapter, we develop 
from Equation (3.18) a system which involves two feedback circuits to account for 
the previous values involved in delta modulation and convolution. We call this the 
reduced delta convolution (RDC) system. 

The RDC system is a particularly simple implementation of Equation (3.18) 
involving only gates and two registers the lengths of which are n\ and n\ + 1 
stages where n\ is the number of bits used to represent the system impulse response 
h(t). The form of the RDC system can be developed by considering the terms of 
Equation (3.18) at times T 3 ,2T 3 and 3T, as follows: 
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where 



y(T.)-C,-[*o-/io]ftr[r,] 
y(2T„) = Cl • [/to • (2x 0 + Xi) + Ai ■ x 0 ] 

= 2r[r,]+r[2T„] 

= y(T,) + I[2T,] 

r [ 2 T 3 \ = C\ • [ho • X\ + h\ • xq] 

J[2T,) = r[T.]+r[2T,] 

= I[T.]+r[2X,] 

/[T,] = r[T.] 

and 

y (3 T 3 ) = Ci • [/io • (3xo + 2xi + £ 2 ) + ^1 * (2xo + ^1 ) + ^2 ■ £ 0 ] 

= 3r[T 3 ] + 2r[2T 3 ] + r[3T 3 ] 

= y( 2T 3 ) + /[3T 3 ] 

where 

r [3T 3 ] = Ci • [/i 0 • x 2 + /ii • xi + h 2 • x 0 ] 

/[3T 3 ] = r[T 3 ]+r[2T 3 ] + r[3T 3 ] 

= /[2T 3 ] + r[3T 3 ] 

Note that the present value of the output is the previous value plus an in- 
crement I. Further, the present value of I is the previous value plus a residual 
r which depends on the present contents of the input register. So, the new data 
updates r which then updates I which then updates the output. These updates are 
accomplished with feedback which involves storage of a discrete voltage as shown in 
Fig. 4.1 where y (j • T 3 ) ~ y(t) , and clock rate = Although Fig. 4.1 has 
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input register 




stores read in on (k-l)th tick and read out on (k)th tick 



Fig. 4.1. Reduced delta convolution system. 
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been developed using transient or initial values (while the input register is becoming 
occupied), it is shown in Appendix B that this same form prevails when the input 
register is filled with values of x(jT 3 ). So, the output of Fig. 4.1 represents the 
convolution of x(t) and h(t ) when the input register is fully occupied, and the 
output is a transient response when the register is not fully occupied (loading and 
unloading). This transient behavior is part of any filter realization that uses delay 
elements (digital or recursive filters). 

A. CIRCUIT REALIZATION OF TRANSVERSAL FILTERS 
USING THE RDC METHOD 

This section first considers a circuit based on Equation (B.ll) which functions 
as a convolver (linear filter) and then a similar circuit which is a correlator. We 
show that a convolver or a correlator can be constructed by using two shift registers, 
ni XNOR gates, one or two amplifiers, one or two summers and two feedback loops. 
We show that by simply reversing the contents of the impulse response register, a 
convolver, in essence, becomes a correlator. 

1. RDC Transversal Filter as a Convolver 

From the preceding discussions and Appendix B, Fig. 4.2 is a diagram 
of a system which can be used to implement a RDC transversal filter which is a 
convolver. 

In the realization of the filter of Fig. 4.2, we show the necessary delays 
as consisting of charge-coupled devices (CCD). Not shown in Fig. 4.2 is the clock 
of rate jr ticks/second used to transport the bit stream through register R x and 
used to transport the multilevel voltages through the CCD. 
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Fig. 4.2. RDC transversal filter convolver. 



62 



In Fig. 4.2, 

conv(j ■ T s ) = output = y(j ■ T s ) ~ y(t) 

Aconv(j ■ T s ) = conv(j ■ T s ) — conv([j — 1] • T s ) 
A2conv(j ■ T s ) A Aconv(j • T s ) — Aconw([; — 1] • T s ). 



In Appendix B, it is shown that 



•n, -1 



ni —1 



A2conu[j • T s ] = Ci • < ^ < Rh(i) > • < R x (i) > - < R x (ni) > • ^ h t 

l i=0 i=0 

where < Rh(i) > represents the contents (±1) of cell i of the impulse response reg- 
ister, etc., and where < R x (i ) > changes with time. So. the value of A2 conv[j • T 1 ,] 
depends on the contents of register R x at time j ■ T s . 

There are several constants in Fig. 4.2 which are defined as follows: 



n, -1 



= amplifier voltage gain = —C\ ■ h x 



i=0 



= amplifier voltage gain = C\ 



A i , Ao , A 3 and A 4 are voltage summers. 

2. RPC Transversal Filter as a Correlator. 

When the RDC system is used as a correlator, then the content of the 
impulse response register is reversed as shown in Chapter 2. Following the same 
procedures given in Appendices A and B, the following result is obtained 



m -l 



n, -i 



A 2 corr(j ■ T s ) = -C\ • < ^ < Rh(i) > ■ < R x (i + !)>-< R x ( 0 ) > • ^ h t 



i—Q 



i = 0 



(4.2) 



where 



A2 corr(j ■ T s ) = Acorr(j ■ T s ) — A corr([j — 1 ] • T, ) 
A corr(j • T s ) A corr(j ■ T s ) — corr([j — 1 ] • T s ) 



63 



and where A'2corr(j -T s ), A corr(j -T 3 ) and corr{j -T s ) are the incremental change in 
the change of the discrete correlation values, the change in the discrete correlation 
values, and the result of correlation at time j ■ T s . 

Fig. 4.3 is a block diagram of a correlator as given in Equation (4.2) 



where 



n 1—1 



— Ci m ^ ^ hi — — fci and 1 4 — — C i — — k'l 



i = 0 



As expected, there is similarity of structure and commonality of hardware 
in the realizations of the convolver and correlator. In fact, if the input register is 
bidirectional and if the amplifiers become inverting, then a convolver becomes a 
correlator and vice versa. 

We can draw some conclusions concerning hardware from Fig. 4.2. 

a. When h(t) begins and ends at 0 volts, then* 

n i—l 



E = 0 



«=0 

and we will not need the amplifier of gain k i nor the summer Ao. 

b. The needed storage equals (2 ni + 1 ) binary stages or cells and 2 single 
stage CCD’s. 

c. The needed number of binary multipliers (XNOR gates) is n\. 

d. The needed number of operational amplifiers is 4. 

e. Ci is a scale factor (could be apart from the IC realization of the RDC 
system). 

f. The same IC can be used as a convolver or a correlator. 

* Given h(ni ■ T s ) = 0. We know h 3 ( 7 Zi • T 3 ) % h(n \ ■ T s ). But 

Tl i — 1 

h 3 (ni ■ T s ) = Ah- hi when h( 0) = 0. So, when h{i%\ ■ T s ) = 0, then 

1=0 

rii — 1 

^2 = 0 , under the condition /^( 0 ) = 0 , (/i a ( 0 ) is always 0 ). 

;=o 
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< R (0) > = h 




j** 



<R (0) > — x 



j +r » i ■ 
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Delay 
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Delay 





output 



A corr( [j+n 1 * 1 ] . T ) one stage CCD | one stage CCD 

corr( [j+n ] - 1] . T ) 



Fig. 4.3. RDC transversal filter correlator. 
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We can now identify some important features of the RDC system. 

a. Uses digital technology (compatible with VLSI). 

b. There axe no code words and hence no synchronization requirements. 

c. There axe no multipliers (XNOR gates accomplish an equivalent function). 

d. Real time operation (no off-line processing). 

e. Easily programmable to represent impulse response of interest. 

f. No limitation on the time duration Td of the signal in matched filter 
applications. 

g. No theoretical limitation on the number of register stages (taps); filter 
sections can be easily cascaded using several IC’s to realize longer 
duration impulse responses. 

h. Reduced hardware requirements (complexity) compared with other 
FIR filter realizations. An example is provided in Chapter 5. 

In Chapter 7, we conclude from simulation results that when used as a matched 
filter or ordinary lineax filter, the RDC system has noise performance which is 
better than predicted by ordinary filter theory. The reasons for this are explored 
in Chapter 7. 

This chapter (4) presents a system (block diagram) which realizes an approx- 
imation to the convolution integral. The system can be used as a linear filter or 
a correlator. The next chapter details the hardware (circuit diagram) required to 
realize the RDC system. 
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V. HARDWARE REALIZATION OF THE RDC SYSTEM 



A. INTRODUCTION 

In this chapter we introduce a possible circuit to implement the RDC system. 
We also compare the needed hardware and complexity of the RDC system and 
a conventional digital matched filter (DMF). We conclude the RDC system has a 
simple layout for IC realization and requires less than 10% of the transistors needed 
for a typical DMF having comparable performance. 

B. CIRCUIT DIAGRAM OF THE RDC SYSTEM 

Fig. 5.1 shows a possible circuit having as output A2 conv(j ■ T s ) which is 
defined in Equation (B.ll) and shown in Fig. 4.2. Register Rh contains the bit 
sequence representing the delta modulation equivalent of the impulse response of 
the filter being realized. Register R x contains (ni +1) bits of the delta modulation 
representation of the input voltage to be filtered. The XNOR gates produce the 
binary products of the contents of register 7?/, and register R x . From Equation 
(B.ll), we need to sum the outputs of the XNOR gates. But, it is difficult to sum 
voltages. Therefore, we use field-effect transistors (FET’s) as voltage- to-current 
converters and then sum currents. Each FET will be in one of two states: 

1) Cut off when its gate has a “—1” (or 0) voltage which means zero current 
passing through it. 

2) Turned on when its gate has a “1” voltage which means I amperes passing 
through it. 
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Fig. 5.1. A circuit which provides A2 conv(j ■ T 3 ). 
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All the currents from the FET’s enter node B to form Itotal leaving node B. 



So, from Fig. 5.1, 



Vi =r 3 - Itotal , 



or, 

V 1 =r 3 -I-(n a ) 

where n a = the number of bits for which <Rh(i ) >=<R x (i ) >; i = 0, 1,2, ... — 1 

Let A = n a — number of bits for which < Rh{i ) >^< R x (i ) >; i = 0, 1, 2, . . . n\ — 1 
So, A = n a — (ni — n a ) = 2 n a — ni\ n a = 0 , 1 , 2 , . . . rt \ . 

As a result of this, 

Ts • /• 

Vi = (A + ni); A = -ni,-ni +2,...,ni -2,r»i (5.1) 

and, 

(Ul )mn — 7*3 • I • Tl\ . 

Assume that ( V\) max = 2V CC . Then V cc = r3 '” 1 — where V cc is an externally 
supplied constant (dc) voltage. Using this in Equation (5.1) gives 



Fig. 5.2 shows V\ as a function of A. 

We can use Equation (5.2) to define V 2 in Fig. 5.1 as follows 



(5.2) 



V 2 = V cc -V 1 = V cc - 



r 3 • I 



A 4- V CI 



r 3 ' I 



Fig. 5.3 shows — V 2 as a function of A. 

/m-i \ 

Let k\ = C\ • I h* ) = constant for the impulse response of interest. From 
Fig. 5.1 



V3 = — V 2 — fcj* < R x (ni ) > 

= -^7— A — ki - < R x (ni ) > 



(5.3) 
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Fig. 5.3. — V 2 as a function of A. 
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From Equation (B.ll) 



A2 conv(j ■ T 3 ) = C i 



{ ni-l 

< R h(i) > ■ < R x(i) > 



- fci- < R x{n\) > 



But, 

Tl\ — 1 

> • < R x (i) >= A 

j=o 

Therefore, we can write 



A 2 cont>( - ; • T s ) = Ci • A — < R x (n\) > • (5-4) 

From Equations (5.3) and (5.4), we can make V 3 = A2conu(j • T s ) if we set Ci = 

laA 

2 * 

Knowing the step sizes Ax and A h, the sampling frequency and the im- 
pulse response register contents and length ni, we can determine the constants C 1 
and k\. Knowing n l5 C\ and choosing a convenient value for /, we can determine 
r 3 and V cc . We choose Vdd in Fig. 5.1 sufficiently larger than ( V 1 ) max so that the 
current I is constant. 

The circuit of Fig. 5.1 has as its output A2conv(j ■ T s ). To find y(j ■ T a ), 
the value of convolution, we have to use two feedback circuits as mentioned in 
Chapter 4. Fig. 5.4 shows one circuit which provides y(j ■ T a ) when its input is 
A 2 conv(_; • T 3 ). The values of r 5 are arbitrary. From Fig. 5.4, 

y(j • T s ) = conv(j • T s ) 

= conv {{j — 1) • T 3 ) + A conv(j ■ T s ) 
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where 



Aconv(j • T a ) = A conv((j — 1)T 3 ) + A2 conv(j • T a ) 
and A2 conv(j • T a ) is given in Equation (5.4). 

C. HARDWARE COMPARISON OF THE RDC SYSTEM AND A 
TYPICAL DIGITAL MATCHED FILTER (DMF) 

We can estimate the hardware requirements of the RDC system using Figures 
5.1 and 5.4. There are four different circuits identified in those figures. There is 
also a delta modulator which could be part of the system or could be off the chip 
when using IC technology. 

Assume 6 transistors are needed to realize each stage of a shift register. Then 
the RDC circuit requires 6ni transistors for the impulse response register and 
6(ni + 1) for' the input register for a total of 6(271! 1) transistors. 

Assume each XNOR gate requires 10 transistors. Then the array of gates 
requires lOni transistors. 

There are ni field-effect transistors. 

Assume 20 transistors are required for each op-amp for a total of 120 transistors. 

The total hardware requirements, neglecting resistors and capacitors, is (23ri! + 
126) transistors plus 2 single stage analog storage elements (CCD’s) and possibly 
a delta modulator. The hardware (in fact, the circuitry) is the same for the re- 
alization of a convolver (lowpass filter, bandpass filter, etc.) as for a correlator 
(matched filter). 
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A conventional way to implement a DMF is presented by Turin [10] (see Ap- 
pendix C). We will compare the hardware required for this DMF with that previ- 
ously tabulated for the RDC system. 

Assume that the sampling frequencies for the conventional method given by 
Turin [10] and the RDC method are 2f max and 2 • F • f max respectively, where 
fmax is the maximum frequency component of the band limited version of the 
input signal and F is the oversampling ratio. As a result of this 




where ni is the register length of the RDC system and L is that of the conventional 
DMF register. 

We can now list the needed hardware to implement a conventional DMF and 
that needed to implement an equivalent RDC filter. The results are in Table 
5.1. 

If we assume each adder needs only 6 transistors, then the needed transistors 
to implement the DMF filter, neglecting the D/A and A/D converters and the 
amplifiers required are 

Q(N + M) • L + 10 M ■ N ■ L + QM ■ N(L - 1) 
where the needed transistors for the RDC system are 

23nx + 126 = 23(F • L) + 126 



Define R? as 



Rj' 



the needed transistors when we use RDC method 
the needed transistors when we use DMF method 
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TABLE 5.1. COMPARISON BETWEEN THE DMF CONVENTIONAL 



METHOD AND THE RDC METHOD 



Conventional DMF 



RDC 



1. N shift registers having Two shift registers having 

N ■ L binary stages and one 2n\ + 1 binary stages 

„ having M • L binary stages 



2. M • N ■ L binary multipliers 
(XNOR gates) 



rii binary multipliers 
(XNOR gates) 



3. M • N ■ (L — 1) summers 



6 operational amplifiers, 
each having less than 20 
transistors. 



4. A/D converter 

5. D / A converter 



Delta modulator 
none required 



6. M • N sets of interconnections n\ interconnections 
where the number of these sets 
equals to L 



For large values of L , then 



A 



R'j' 

common value for M and N is 



23F 

6(iV + M) + 16 M ■ N 

8 



(5.5) 



then Rt « — (5.6) 

1 1120 V ; 

From Equation (5.6) we conclude that the circuit complexity is equal for F « 50. 

When L = 1000, M = N = 8, and F = 4 then Rt = 0.0823. In general, we 
conclude that the number of transistors needed to implement the RDC method 
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is less than 10% of that needed to implement the conventional DMF filter. In 
addition the circuit layout of a conventional DMF method is more complex than 
that of the RDC system. 
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VI. PERFORMANCE OF THE RDC METHOD 



In this chapter we present simulation results which show the transfer function 
of a RDC lowpass filter (LPF), the response of a RDC LPF to two signals, one in 
the pass band and the other outside the pass band, and the output of a RDC LPF 
with a square wave input. Included are responses of filters matched to a raised 
cosine pulse and to a chirp signal. 

A. LOWPASS FILTER PERFORMANCE 
1. Transfer Function of the RDC System 

Fig. 6.1 is a block diagram of the RDC system used to find the transfer 
function when used as a LPF. We assume a LPF cutoff frequency f c of 1.0 Hz. 
The impulse response h(t) is, then, 

sin(27 rf c t) _ 0 sin(27rt) 

' /c ' 27 rf c t 27 Tt 

where the corresponding causal impulse response h c (t) together with its staircase 

representation h cs {t ) are plotted in Fig. 6.2. 

We obtain the transfer function by applying to the RDC system the signal 

x{t) = Ai • sin(27r fit). 

Following the transient time, which is dependent upon the register length 
n\ and the sampling interval T,, the output y(t ) is 

y{t) = Aq • sin(27r/;< + <f>) where <f> « 0. 

The ratio ^ as a function of frequency is taken as the magnitude of the RDC LPF 
transfer function. 




77 



analog 

input 

x(t) = A 



/ 



delta modulator 




RDC 


output 




LPF 


y(t) 


t ) 


H(f) J 


{ 





- f 



Fig. 6.1. Block diagram of the system 
used to find the RDC LPF transfer function. 
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6.2. The causal transfer function for the LPF h c (t) 
together with its staircase representation h C3 (t). 



The register length n x depends on T s and Td which is shown in Fig. 6.2. In 
fact, rii = For a fixed T,, increasing n\ permits a more inclusive representation 
of the impulse response h(t). In the following, we measure Td in terms of the number 
of zeroes about the main lobe of h c (t ) and included in Td . 

Figures 6.3, 6.4, and 6.5 show the magnitude of the transfer function for 
the RDC LPF where the number of zeroes of the impulse response are 2, 6, and 8 
respectively. As expected, the LPF transfer function “improves” as more zeroes of 
h c (t) are included. 

Fig. 6.6 shows the magnitude of the transfer function of the RDC LPF 
when Td includes 6 zeroes of h c (t ) and when rectangular and Hamming windows 
are used. We see that the ripple in the transfer function is reduced using the 
Hamming window. 

When the step size Ah becomes larger than the side lobes of the impulse 
response, then the delta modulator considers the side lobes are of equal amplitudes 
(see Fig. 6.2). To alleviate this problem we can choose a smaller step size and 
increase the sampling rate to maintain the DM slope-following capacity. This 
results in increased register length n^. An alternative is use of a variable step size 
RDC system which is treated in Chapter 7. 
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Fig. 6.3. The RDC LPF transfer function when Tj includes 2 zeroes of h c (t). 
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Fig. 6.4. The RDC LPF transfer function when Tj includes 6 zeroes of h c (t). 



Mh CQ 
'w' TU 




83 



Fig. 6.5. The RDC LPF transfer function when Tj includes 8 zeroes of h c (t). 
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Fig. 6.6. The RDC LPF transfer function when Tj includes 6 zeroes of h c (t) and 
when rectangular and Hamming windows are used. 



2. The Output of the RPC LPF When the Input is the Sum 
of Two Sinusoids, One in the Pass Band and the Other 
Outside the Pass Band 

Assume that the RDC LPF with f c = 1.0 Hz, has an input 



x(t ) = 0.5 sin 




+ 1.0 sin(u;it) 



where izq = 2irfi and f\ = 1.2Hz. 



Figures 6.7 and 6.8 show the input x(f), its staircase representation x s (t), and the 
output y{t) of the filter. We see that the 0.4 Hz sinusoid is preserved in the output 
while the 1.2 Hz term is suppressed. The RDC system when configured as a LPF 
does discriminate (filter) on the basis of frequency. 

3. The Output of a RDC LPF with a Square Wave Input 

Assume f c = 1.0 Hz. Assume the input x{t) is a square wave of frequency 
/o = 0.25 Hz, where Fig. 6.9 shows the input x(t) and its staircase representation 
x s (t). It is expected that the first and third harmonics of the input will appear in 
the output of the LPF and all other harmonics will be suppressed. Further, the 
amplitude of the third harmonic should be one third that of the first harmonic 
according to the Fourier series representation of the square wave. The simulation 
result agrees with this expectation as shown in Fig. 6.10. 

Assume that the square wave input has a fundamental frequency of 0.5 
Hz. Fig. 6.11 shows that all harmonics are suppressed in the output of the RDC 
LPF, as desired. 
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Fig. 6.7. The input x{t) and its staircase representation x a (t)\ 
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Fig. 6.8. The output of the RDC LPF when the input is the sum of two sinusoids, 
one in the pass hand and the other outside the pass band. , 
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Fig. 6.9. Input square wave x(t) and its staircase representation x 3 (t). 
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Pig. 6.10. The output of the RDC LPP when the input is a square wave 
having a fundamental frequency of 0.25 Hz. / 



o 
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Fig. 6.11. The output of the RDC LPF when the input is a square wave 
having a fundamental frequency of 0.5 Hz. / 



B. MATCHED FILTER PERFORMANCE 



1. Filter Matched to a Raised Cosine Pulse 

Fig. 6.12 shows a block diagram of the system used to simulate the im- 
plementation of a filter matched to a raised cosine signal. The impulse response of 
the matched filter 



h(t) = x (T d - t) = 



A (l +cos (2 tt - |))) 0<t<?f 

-A (l + cos (2w ^ <t <T d 



where a:(t) is the input signal and Tj, is the time duration. Fig. 6.13 shows the 
raised cosine impulse response where A was chosen to be equal to 1.0 volt and 
Tj, = 1.0 second. 

Figures 6.14 and 6.15 show the results of correlation when we use the 
conventional methods [y (t)] and the RDC method [y(t)] with an oversampling 
ratio F = 8.0. Note that these results have the same shape. The theoretical 
maximum of the correlation is 



\a 2 T a = 1.5 

Simulation results give a maximum value of correlation equal to 1.49999 when using 
the conventional method and 1.54536 using the RDC method. 
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Fig. 6.12. Block diagram of the RDC matched filter that was simulated. 
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Fig. 6.13. Raised cosine MF impulse response. 
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Fig. 6.14. Output of the raised cosine MF when conventional methods are employed. 
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Fig. 6.15. Output of the raised cosine MF when the RDC method 

is employed when F = 8. i 



2. Chirp Matched Filter 



The chirp or linearly frequency modulated signal has been popular in pulse 
compression applications in the past because it is possible to realize filters matched 
to the chirp signal using acoustic devices. The instantaneous frequency fi(t ) as a 
function of time is given in Fig. 6.16 where /,•(<) = f\ + 1 ^ t where T<* is the 

chirp pulse time duration. So, the chirp signal (voltage) is then equal 



A • sin 



(^Wit, + 



(W2 ~V>l) 2 \ 

2 T d ) 



A block diagram of the system simulated is shown as Fig. 6.17. Figures 
6.18, 6.19, and 6.20 show the input signal x(t), the output y cvi (t) when we use 
a conventional method, and the output y(t) of the RDC matched filter. System 
parameters used are fi = 0.0 Hz, /2 = 5.0 Hz, Tj = 4.0 seconds, A = 1.0, F = 4, 
and n\ = 160. 

From Figures 6.19 and 6.20 we observe that the maximum of y cvi (t ) is 
1.88820 and of y(t ) is 1.84669 where the theoretical maximum is 1.8879557. After 
the input expires the output is constant (non-zero). This dc offset is explained 
in Appendix D. From Equation (D.3) the dc offset is expected to be equal to 
±4p (total area under h(t)). 

Increasing the oversampling ratio F to 32 decreases the dc offset. See 
Figures 6.20 and 6.21. This is so because when F increases, Ax decreases to keep 
the value of the slope-following capacity constant. And when Ax decreases, then 
from Equation (D.3), the dc offset decreases. 

Note that in the case of the raised cosine matched filter, the result of 
correlation has no dc offset because the total area under h(t) equals 0 (see Fig. 
6.15). 
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Fig. 6.16. The instantaneous frequency of the chirp signal. 
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Fig. 6.17. A block diagram of the RDC chirp matched filter that was simulated. 
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Fig. 6.18. The input chirp signal. 
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Fig. 6.19. The output of the chirp MF when using conventional methods. 
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Fig. 6.20. The output of the chirp MF when using the RDC method for F 
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Fig. 6.21. The output of the chirp MF when using the RDC method for F = 32. 



VII. NOISE PERFORMANCE OF THE RDC SYSTEM 



A. INTRODUCTION 

A DM converts an analog signal to a binary signal. This implies quantization. 
As a result of this and because of noise which may appear on the input, the output 
is subject to distortion. We identify two types of distortion and call their causes 
self noise and external noise. Self noise is caused by the delta modulation process. 
External noise is present at the input to the DM. 

In this chapter we consider the effects on the output of the RDC system of 
first self noise and then external noise. Simulation results for lowpass filters and 
matched filters are presented. The idea of filtering on the basis of the slope of a 
voltage is introduced. Also included is the variable step size RDC method used 
to overcome a self noise effect on the representation of the impulse response of a 
filter. 

Self noise is affected by the step size and sampling rate used in DM. Therefore, 
before considering self noise, we first establish an optimum step size. 

B. OPTIMUM STEP SIZE 

In this section we define and derive the optimum step size for use with the 
RDC method. 

Let x(t) be a zero-mean random input to the two level quantizer shown in Fig. 
7.1. (In Fig. 3.1, the hard limiter, amplifier of gain G , and the flip flop can be 
replaced by a quantizer in the sense given in Fig. 7.1). 
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Fig. 7.1. Two level quantizer, (a) Block diagram, (b) Transfer characteristics. 

(c) Example of input pdf. 
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Let x(t ) have a probability density function (pdf) p x (-) and variance o\. (If 
the input x(t ) has a nonzero mean, we can subtract it from the input and add it 
back after quantization.). The quantizer output y is 

+Ax if x G (0, oo) 

(7.1) 

—Ax if x € (— 00 , 0) 

The quantization error q e = x — y is a random variable with pdf p q (-\ and 
variance cr 2 , where 



y = Q( x ) = 



But, 



So, 



OO 

Vq = E [ql\ = j ql- Pq(q) ■ dq- 



q t = x — Q(x) is a function of x. 



= E [tf] = J' 



(x - Q(x)) 2 -Px(x) • dx 



/ 

K> 

C 

/ 






x 2 • p x (x) - dx + 2Ax • J x • p x (x) • dx 



+ (Ax) 2 • J Px(x) • dx — 2Ax • j x ■ p x {x) ■ dx 



/ 



If p x (x) is symmetric about zero, then 



a 2 = a 2 — 4Ax 



q “i 



7 ~ x - 

n 



p*(x) • dx + (Ax) 2 



(7.2) 



We define the optimum step size Ax opt as that which minimizes <r 2 . To find Ax op *, 
form 



da 2 

q 



d( A 






x ■ p x (x ) • dx + 2Ax = 0 
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and solve for 



AXopt — 

Using this value in Equation (7.2) gives 



opt — 2 f x j)j^(x) • dx . 



(7.3) 



a 2 = cr 2 x - 2Ax • Ax opt + (Ax) 2 when Ax ^ Ax opt 



(7.4) 



and 



min {a 2 } = a\ - (Ax opt ) 2 when Ax = Ax opt 
If we define a quantizer performance factor as 



(7.5) 



_ G 1 
q ° 2 



(7.6) 



then, 



min 



= mm{^} = _ (Ax opt ) 2 

-r2 1 ^2 






(7.7) 



where the inverse of € 2 can be considered as a signal to noise ratio. 

Table 7.1 lists for different probability density functions and indicates 

the corresponding minimum value of € 2 . When the input is a sine wave, then the 



pdf is an arcsin function and ^° Bt = 0.9 (Jayant [31]). In all simulations this 
value of optimum step size was calculated for the applied voltage and then used to 
obtain the results. 

The value of Ax opt is that which minimizes the variance of the difference be- 
tween the input and the output of the two-level quantizer. The two-level quantizer 
is part of a delta modulator. In fact the DM of Fig. 3.1 is formed by an integrator 
or accumulator in a feedback path from the output of a two-level quantizer to its 
input. We now have a new variable which is the gain g in the feedback loop. In 
the next section, we derive the optimum value of g. 
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TABLE 7.1. OPTIMUM STEP SIZE AND MINIMUM VALUE 
OF THE QUANTIZER PERFORMANCE FACTOR FOR 
DIFFERENT INPUT PDF’S. 



pdf 




mzn {£?} 


Uniform 


v/3 

9 


0.25 


Gaussian 


v/i 


0.363 


Laplace 


l 

72 


0.5 


Gamma 


1 


0.667 



1. Optimum Gain of DM Feedback Loop 

In Fig. 7.2, assume that the accumulator has an inherent 1 bit delay; then, 



x a (n -T a ) = g • x([?r - 1] • T a ) . 



As a result of this 

-x^(n • T s ) = x(n ■ T s ) - x 3 (n ■ T,) = x(n ■ T s ) - g ■ x ([n - l] • T s ) 
(x A (n • T a )f = (x(n • T a ) - g • x{[n -T] • T s )f = 



Therefore. 



= + 9 



~ -9 • Pi ■ 



(7.8) 



where a\ is the variance of the input signal to the quantizer (hard limiter) and p\ 
is the value of the variance normalized autocorrelation function U„(r) of x(t) for 
r = T a . 

Q g.2 

Setting -Qp = 0 yields the optimum value of y and min{c } the minimum error 
variance. 



doj 

dg 



= (2<7 - ) • a\ 



= 0 
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Fig. 7.2. Block diagram of the delta modulator and demodulator. 
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(7.9) 



Therefore, 



ffopt = pi and min{a\j = (1 - p\)crl 



We see from Equation 7.9 that the error variance is less than the input signal 

2 

variance a\ for all non-zero values of pi. We call a prediction gain Gp. 

** A 
2 

The maximum value of Gp = max{Gp} = — ; — - x -- 9 — = (1 — p ?) -1 (7.10) 

Processors of speech and images are often based on long-term autocorre- 
lation functions. In general, for speech pi has a value close to 0.9 (Jayant [31]). In 
sub-optimal design, pi is taken as 1 and so g = 1 . 

The preceding consideration of a two-level quantizer and feedback gain 
can now be applied to determine the optimum step size for DM analog-to-binary 
conversion. 

2. Optimum Step Size for Delta Modulation 

From Fig. 7.2 we can define for a DM signal the reconstruction error 
r(n • T 3 ) as follows 

__ r{n-T 3 ) = x{n-T 3 )-y{n-T 3 ) 



where y(n • T 3 ) is the output of the delta demodulator which is equal to x 3 (n • T s ). 
We can prove that the reconstruction error has a minimum variance 



• r 2 i • r 2 i ™in{e 2 } -(1 - p\) . 
min \cr \ = min \ a a \ = 7 ^ : — r— rr-a 

1 ' ' (Wl-minjsU) 



(7.11) 



where we note that the minimization in Equation (7.11) is with respect to Ax and 
that this result holds for g « pi (Jayant [31]). 

From Equation (7.7) we can say that 



Ax opt — 




min 




= 0 • a 



A 



(7.12) 
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where © is a constant which depends on the pdf of the input x(t). Also, from 
Equations 7.7 and 7.11 



= 



(1 — Pi) 



(1 — p\ • min { €g } ) 
So, from Equations (7.12) and (7.13) 






(7.13) 



Ax 0 pt — 0 • 



(!-/>?) 



(7-14) 



(1 - p\ - min {e\}) 

where mm{6,} < 1, and can be found using Equations (7.3) and (7.7) or Table 
7.1 for the pdf of interest. 

If p\ — ♦ 1 (which is the case for speech and video signals for sampling 
frequencies higher than the Nyquist frequency), Equation (7.14) can be rewritten 
as follows 

^^X opt — * 



(1 -Pi) 



(1 -mm{6 2 }) 



< 7 ,. 



Then using Equation (7.12), 



A x opt = yi - 



Pi °x- 



(7.15) 



For pdf’s of many common input signals, 



max{Gp} = r oc F 

1 — Pi 



Pi 

(Jayant [31]) where F is the DM oversampling ratio defined in Section V.C. 
Let 



max{G P } = Cp g F 2 = 



1 



i -pi 



(7.16) 



where C pg is a constant. From Equations (7.15) and (7.16) 



A x 0 pt — 



Cpg-F 



(7.17) 
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For example, for a sinusoid of frequency f\ and amplitude A which is sampled at 



frequency /, 



x(n • T s ) = A • sin(uq • n • T s ) 



Therefore 



R XX {T 3 ) = A 2 sin(wi • n • T s ) • sin(uq • [n + 1] • T s ) 

^2 . 

= — cos(wi -T s )- cos(wi(2 n + 1) • T 3 ) 

A 2 

= — cos(u; 1 • T s ) 




(7.18) 



and 




(7.19) 



So, 




But 



COS X — \ 1 — - — 

2 4! 



Therefore 




where F » 1. 



So 




As a result of this 




From Equation 7.16 



1 



(7.20) 
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From Equations (7.17) and (7.20) 



^%opt 



7 r 




(7.21) 



But 



Therefore, 







\/2 7T • A ■ fi 

fs 



4.44 A ■ h 

fs 



2.22 A 

F 



(7.22) 



Having established optimum values of step size and gain for a DM, we now 
can use these results to determine values of self noise in the ADC process using 
DM. 



C. SELF NOISE 

Chapter Three identifies two types of self noise: granular (quantization) noise 
and slope overload noise. 

1. Granular Noise Power 

If |x(t)| is less than the step size as shown in Fig. 7.3, the granular noise 
is given by 

e g = x(t) - x 3 (t ) 

Generally, we have granular noise for all x(t) because of quantization. To 
precisely define t g , let us assume that the analog input to the DM modulator x(t) 
and the output of the DM receiver are given in Fig. 7.4. Assume also that the 
step size Ax and the sampling interval T s are infinitesimal, but with their ratio the 
same as the actual system. The output of the DM receiver will then be the smooth 
function x sa (t) shown in Fig. 7.5. 



Ill 



We define slope overload noise e s0 as the difference between x(t) and x s3 (t); 
the remaining distortion ( x ss (t ) — 15(f)), is granular noise. 

Granular noise e g is essentially uncorrelated with the signal (Peebles [32] 
and Protonotarios [33]) and is approximately uniformly distributed (Goodman [29] 
and Protonotarios [33]), over the interval (—Ax, Ax). See Fig. 7.6. 

Therefore, 

E(tg) = J ( 2 ^) • de 9 = °> 

-Ax 

and the granular noise variance is given by 




(Ax) 2 
3 * 



(7.23) 



So, the granular noise power increases as the square of the step size. 

2. Slope Overload Noise 

In general the slope overload noise is correlated with the message (Peebles 
[32]). A number of analyses to find the slope overload noise power are introduced 
by Protonotarios [33], Zetterberg [34], Rice (with O’Neal) [35], and Abate [36]. All 
these results do not either individually or collectively pertain to all slope-following 
capacities and input spectra. Greenstein [30] found an expression for slope overload 
noise power that is accurate for all slope-following capacities and input spectra of 
possible interest for linear delta modulators having Gaussian random inputs. 

Let us define the slope overload factor 




(7.24) 



where x' a is the slope-following capacity. A highly accurate approximation to the 
slope overload noise power <r 2 0 for linear delta modulators having Gaussian inputs 
and for all spectra, that was given by Greenstein [30], is as follows 
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Fig. 7.3. Illustration of the case where |x(f)| < Ax. 
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Fig. 7.4 The input signal and its staircase representation. 




Fig. 7.5. The input signal and the smooth staircase representation. 
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Fig. 7.6. Probability density function of granular noise. 



115 



<j 2 so = cr 2 [ 1 + 2.753s + 2.952s 2 ] exp(-0.341s 2 ) 

• exp {(ax — 2.753)s + a 2 [exp(a 3 s + a 4 s 2 ) — l] } ; 0 < s < 4.0 (7.25) 

where the expression given by Equation (7.25) predicts cr 2 0 with an accuracy of 
1 dB or better for s between 4.0 and 6.5. For s > 6.5, <J 2 so is at least 119 dB below 
<7 2 . The variables a 4 , a 2 , a 3 , and a 4 are determined by the spectrum of the input 
(Greenstein [30]). For example, for bandlimited white noise 



a\ = —0.036, a 2 = 0.37, a 3 = —3.83, and a 4 = —5.9. 



(See Fig. 7.7.) From this section and section VII. C.l we can draw an error variance 
curve for DM (See Fig. 7.8.), where 

= oio +-*; (7.26) 



Note that as Ax falls below the value Ax opt , the distortion increases more rapidly 
than when Ax takes on values greater than A x op (. We say that the delta modulator 
“fits” the signal if Ax = Ax op *. 

3. Signal to Noise Ratio Improvement (SNRI1 Due to Delta 
Modulation 

Assume that we have chosen the step size to be equal to Ax opt where we 
have as input of the RDC system the signal alone. Then, 






( Ax op ( ) 

3 



(7.27) 



where <r 2 x , a 2 SOx , and a 2 x are the total quantization noise power, the slope overload 
noise power and the granular noise power respectively. 
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Fig. 7.7. Normalized slope overload noise as a function of,s. 
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Fig. 7.8. Error variance curve for DM. 



If we have as input noise alone (which has power a\) and the step size is 
A x opt then 



2 2 



+ <7 



3n 



(7.28) 



where > a ^o n •> and a2 g n 3X6 the total quantization noise power, the slope overload 
noise power, and the granular noise power respectively. 

Neglecting the granular noise, which is practically very small compared to 
the slope overload noise, the effective noise to the RDC system will be decreased 
by a \o n • That is > 



the effective input noise power = a 2 n — a 2 SOn = < 7^(1 — *4) (7.29) 



where 

2 

A = < 1. (7.30) 

See Fig. 7.9. From Equations (7.25), (7.29), and (7.30) we can find the input 
effective noise power when we have noise as input. 

When we have as input to the RDC system the signal as well as noise, the 
root mean square value of the slope of the input will increase. From Equation (7.24) 
as ( 37 ) increases (where Ax and T s are not changed), s decreases. From Fig. 
7.7 as s decreases, the slope overload noise power, in this case <J 2 0x+n , increases. 
Neglecting the granular noise, the output power of the RDC system (neglecting 
any SNRI due to matched filtering) is 



Therefore, 



<r 2 x + cr 



2 

n 



— ( 7 



2 

5 °r + n 



(SNR ) 0 
( SNR)i 



<j 2 — <7 2 

u so x + n 




1 - D 



(7.31) 
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Fig. 7.9. Error variance when the input does not fit the delta modulator. 



where 



2 

D = *“»+" ■ < 1 

So, the signal to noise ratio improvement due to using delta modulation = > 1. 

D. THE ERROR IN CONVOLUTION WHEN THE OUTPUT BITS 
OF THE DELTA MODULATOR ARE IN ERROR 

The DM is an ADC. The bit stream output can be transmitted as a baseband 

signal (twisted pair, coaxial cable, or fiber) or as a bandpass signal (modulated 

carrier). Errors may occur in transmission. 

A known signal at the input to a DM creates a known bit stream. Noise added 

to the signal creates a different bit stream. The difference can be considered as bits 

in error. The relationship between the noise characteristics and the bits in error is 

not generally known. 

In this section, we consider the effect on the output of the RDC system of bits 
in error. Because there are two feedback circuits in the RDC system, bits in error 
on the input affect future values of the output. 

From first principles of delta convolution (Appendix A), the error in the con- 
volution due to one bit in error is given by 



= ±2Ci | ^2 hi • ( n i - 0} 



= ±€l 



(7.32) 



where ei is a constant for a given h(t) and T s . 

Assume that we have two bits in error. Then the bits in error can be as follows. 
Both can be +1; both can be —1; or they can be of opposite sign. There are 



~ ) combinations of two l's 

2 



combinations of a single 1 and a single — 1, 



and 



0 



combinations of two — l's. 



m 



where I ) is the binomial coefficient of m things taken n at a time. 



n 



Therefore the error in convolution will be 



H 1 

2ei with probability = - 

f 2 ) 2 

0 with probability -yy- = - 



and 



( ) 1 

—2t\ with probability -^y- = 

Generally, if we have m bits in error and m is even, the error in convolution 
is given by 



m ■ 


with probability 


(SI 

2 m 


(m — 2) • t\ 


with probability 


LSil 

2 m 


(m — 4) • ei 


with probability 


iZA 

2 m 


■ 




2 m 


2e x 


with probability 


0 


with probability 


(m^V) 
2 m 


(m — k)e i 


with probability 


2 m 


and — m • ei 


with probability 


(11 

2 m 
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where k is even number. See Fig. 7.10. 

So, 

Cconv — 0 

and 

a “fn 2 / \ 

™ * } ( 7 . 33 ) 

Jt=o '■ 2 ' 

where m and fc are even numbers. 

In the same way, if we have m bits in error, where m is an odd number, the 
error in convolution is given by 

( m ™i) 

(m — k) ■ ei with probability — ; k = 0, 2, 4, ... , 2m 
So, 

tconv = 0 

and 

4^=|rE(( m - fc ) ie i) 2 (m-f) (7 ‘ 34) 

where m is odd and is even (see Fig. 7.11.). Note that Equations (7.33) and 
(7.34) are exactly the same except that m is even in Equation (7.33), and it is 
odd in Equation (7.34). Note also that the error in convolution, for a certain error 
pattern (where all bits in error have left the first ri\ stages of i? r ), is independent 
of time. 

When the bits in error propagate through the first n\ stages of register R x , 
the error in convolution is dependent upon the position of the bits in error and so 
it is dependent on time. We want to find the signal to noise ratio at those times 
when the bits in error have passed through register R x \ . Assume that we have m 
bits in error. Therefore, from Appendix A, 



123 




Fig. 7.10. The probability of error in convolution, m even. 
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Fig. 7.11. The probability of error in convolution, m odd. 
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the output « C\ 



m -l 



y: hi(ni - i ) 



i=0 



jf-T»l -1 



I] x ' 



i=o 



+ noise at t = j • T a 



where the noise has zero mean and variance given by Equations (7.33) and (7.34). 
Therefore, 



SNR 



j-n i -1 

^ E 

i=0 



2c 



2 m — 2 






k=0 



j-n t-1 

E X i 
1=0 



m— 2 






fc=0 



(7.35) 



Assume that the probability of a bit in error is p e . So, the average time period 
between two bits in error equals ^ L . In practice p € < 10 -3 . As a result of this the 
expected average time between two bits in error is > 1000T S . 

From previous discussions, we know that the effect of a bit that has passed 
through register R x i is the same as adding dc voltage to the output. As a result of 
this and because the average time period between two bits in error is greater than 
1000T,, an ac coupling in the output of the RDC system will tend to cancel the 
effect in the output of the bits in error. 



E. SIMULATION RESULTS 

In this section we present simulation results when the RDC system is used to 
implement a chirp matched filter and a linear filter (lowpass filter) when additive 
noise is part of the input. 
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1. RPC Chirp Matched Filter 



The system simulated is shown in block diagram form in Fig. 7.12. The 
input noise to the lowpass filter is white noise with zero mean and standard de- 
viation sd\. The bandwidth of the LPF equals that of the chirp signal, and the 
output band limited noise of the LPF has a standard deviation sd. 2 , where 

x(t) = A sin ) . C-7-36) 

f 2 = 5.0Hz, Td = 4.0 seconds, .F = 8, and .4 = 1.0 

When the band limited noise has a standard deviation sd ,2 equal to 0.75, the average 
results of six computer runs (where each run had a different seed for the random 
number generator) are as follows. The maximum value of correlation = 1.641; the 
output noise has approximately zero mean and a standard deviation 0.2768 where 
the average chirp signal power numerically is 0.472. See Figures 7.13 and 7.14 for 
the input and the result of one run. From the given results 

0 472 

(SNR) , = pjjj = -0.762 dB 
< SN «>° = = 15A59dB 

Therefore, the signal to noise ratio improvement SNRI = 15.459 — (—0.762) = 
16.221 dB. 

Fig. 7.15 shows the SNRI for different values of sd 2 where theoretically 
the SNRI = Td • B, and where Td • B is the pulse duration-bandwidth product 
of the chirp pulse. (See Horrigan [11].) Hence, for this example, the theoretical 
SNRI = 4(5) ==> 13dB. So, SNRI using the RDC method exceeds the maximum 
theoretical value by about 3 dB for this case. 

Denyer [37] introduced a digital technique called Deltic processing to im- 
plement chirp matched filters. We found that the SNRI using the Deltic processing 



127 




band limited noise 



Fig. 7.12. Block diagram of the RDC chirp matched filter that was simulated. 
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Fig. 7.13. The input of the RDC chirp matched filter. 
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Fig. 7.14. The output of the RDC chixp matched filter. 
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Fig. 7.15. SN RI as a function of the standard deviation of the band limited noise 

when using the RDC system. 



technique to implement a chirp MF is less than that obtained by the RDC method 
by 9 dB. 

The noise performance of the RDC system can be explained as follows. In 
fact, there exists a double match. The first is due to the matching of the filter to 
the signal. The second is because the delta modulator step size is adjusted to fit 
the input signal so there is no slope overload. When noise is added to the signal, 
slope overload occurs and the delta modulator does not always track the signal 
plus the noise. This lessens the effect of the noise input. 

2. Lowpass Filter 

In Fig. 7.16, assume that the LPFs have a cut-off frequency f c equal to 
5Hz and that fi equals 2.0 Hz where the input is 

i(t) = x{t) -f n(t) 

where n(t) is zero mean white Gaussian noise with variance = (0.5) 2 . Assume also 
that only 5 zeroes at each side of the main lobe of the impulse response of the 
filters are considered. 

Figures 7.17, 7.18, 7.19, and 7.20 show the input i(t), and the outputs 
j/i (t), J/2 (^) > and yz{t). Fig. 7.18 shows that the DM acts like a filter because the 
step size for the delta modulation process is adjusted to fit the signal. The slope 
factor is s = 1.0. When noise plus signal are applied to the input of the delta 
modulator, the slope factor s decreases leading to slope overload distortion. This 
actually reduces the noise effect. In fact, the effective noise power into the RDC 
system is less than the input power as derived in Equation (7.29). 

We can compare the noise performance of the RDC filter with a conven- 
tional LPF by defining 
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Fig. 7.16. Block diagram of the system used to simulate the outputs of 
the delta demodulator, the RDC LPF, and the 
conventional LPF. 
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Fig. 7.17. The input to the RDC and conventional lowpass filters. 
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7*18. The output of the delta demodulator. 
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• 7.19. The output of the RDC lowpass filter. 
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tig. 7.20. 7 lie output of the conventional lowpass filter.' 



d 2 {t) = yi{T tr + t) - x(t ) 



and 

dz{t) = y 3 (i) - x(i) 

where T( r ~ 1 second is the transient time needed to fill the n\ stages of the register 
R x i of the RDC system. 

The simulation results show that the average power of d 2 (t) is equal to 
0.01324 where that of d 3 (i) is equal to 0.01902 which means that the output of the 
RDC LPF is closer to the original signal than the output of the conventional LPF. 
We then conclude that the “effective” noise into the RDC system is less than the 
input noise of the conventional LPF. 

Another way to show that the delta modulator is a “slope” filter is to use 
the block diagram given in Fig. 7.16 to find ri,r 2 , and r 3 versus a\ where 




and where aj, a^cr^, and erf are the variances of yi(t),y 2 (t),y 2 (t) and i(t) respec- 
tively. We consider the case where the input i(t) is white Gaussian noise only and 
where the step size of the delta modulator is chosen to fit a sinusoid of amplitude 
1 volt and frequency equal to the bandwidth of the LPF which is assumed to be 5 
Hz. 

Simulation results are summarized in Fig. 7.21 as plots of ri,r 2 , and r 3 
versus cr?. Note 
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Fig. 7.21. The values of ri,r 2 , and r 3 as a function of a\. 



• For the RDC LPF, rj and decrease as the variance of the input increases. 
This is because an increase in the variance of the input creates additional slope 
overload conditions. The delta modulator then acts as a “slope” filter. The 
slope of the output of the delta demodulator increases but at a smaller rate. 
As a result of this the values of r\ and 7*2 decrease as a\ increases. 

• When a\ is small compared to the step size, granular noise occurs and so 7*1 
will be greater than 7-3. 

The simulation results verify that the noise performance of the RDC sys- 
tem is better than that expected from the theory of linear filters. The reason for 
this is that the DM portion of the RDC system filters on the basis of slope as well 
as frequency of the applied voltage. 

F. USE OF VARIABLE STEP SIZE WITH THE RDC METHOD TO 
OVERCOME GRANULAR NOISE OF THE IMPULSE 
RESPONSE REPRESENTATION 

Assume that we want to represent the impulse response h{t) for a LPF or a 
band suppression filter using DM where the step size A h is constant. For a given 
sampling rate, if A h is appropriate for the main lobe (to prevent slope overload) 
then the DM values for all side lobes whose amplitudes are less than the step size 
are successive values of +1 and —1. This means that the side lobes are represented 
as being constant. On the other hand, if we choose the step size A h appropriate 
for the side lobes, the DM cannot track the main lobe without increasing the 
sampling frequency. This requires an increase in the number of stages of the two 
shift registers Rh and R x . One solution is to use different step sizes for different 
regions of h(t). We call this the variable step size RDC method. 
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1. Implementation of the Variable Step Size RPC Method 

To implement the variable step size RDC method, we first define h(t) in 
the range seconds. To make the impulse response causal, shift the time axis 

by seconds. An example is shown in Fig. 7.22. 

Now divide the causal impulse response into M different time intervals so 
that in each interval a single step size is appropriate. In Fig. 7.22 M — 5 and, 
because of symmetry, the step size for region A equals that of E and the step size 
for region B equals that of D. Assume that the M intervals have ni, n 2 , . . . , 
samples and step sizes respectively. 

Using the procedure given in Appendices A and B, it is easy to show that 

{ 7M-1 

y: [< Rh(i ) > • < Rx(i) >] • A hi 

i= o 

ni+n 2 -l 

+ ^2 [< Rh(i) > • < Rx(i) >] ' A/l 2 

i=n l 

n i -hn2 -hri3 — 1 

+ 53 [< ^*(0 > • < ^ r (0 >] • A h z 

i=n 1 -f-n 2 

+ 



+ 



iV-l 



+ ^ [< Rh(i ) > • < R x (i ) >] • A/iAf 



i=N-n 



M 



where 



and 



- < R X (N ) > -k 5 >A x-T s 



M 






Hi -1 



2=1 



n l + n 2 ”1 



N - 1 



(7.37) 



ka — 



A/ii • y ^ hi + a/i 2 • y ' /it + + a/im ■ y ^ hi 

«=0 i=n i i= N — n m 

constant for a given impulse response and given step sizes. 
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Fig. 7.22. The causal impulse response of a band suppression filter 
divided into 5 regions, each of which has an appropriate step size. 



Equation (7.37) can be implemented as shown in Fig. 7.23 where the constant 
Ax • T s can be off the IC chip. 

2. Applications 

This section presents simulation results obtained using the variable step 
size RDC method to implement a lowpass filter and a band suppression (notch) 
filter. 

a. Variable Step Size RDC Lowpass Filter - 

Assume that the cutoff frequency of the LPF is 1.0 Hz. Figures 7.24. 
7.25 and 7.26 show the impulse response and the staircase representations ( h v (t ) 
and h c t(t )) when we use variable and constant step sizes and where we take into 
consideration 40 zeroes on each side of the main lobe. The number of intervals is 
M = 21. Fig. 7.27 shows the resulting transfer function when we use variable step 
size. Also shown is the transfer function when a constant step size is used. Using 
the variable step size RDC method, the ripple in the transfer function is reduced 
and the filter skirts axe steeper. 

b. Variable Step Size RDC Band Suppression (Notch) Filter 

Fig. 7.28 shows the transfer function of a band suppression filter where 

f b = 1.0Hz, / a = ^ = 0.5, and Af = 0.375. 

<u 

and where 

H„(f) = H,(f ) - H„U) 

Figures 7.29, 7.30 and 7.31 show the impulse response and the corresponding stair- 
case representation when we use the variable step size and constant step size RDC 
method. The number of intervals for the variable step size is M = 17. The transfer 
functions obtained are shown in Fig. 7.32. 

From Figures 7.27 and 7.32, we can say that the variable step size 
RDC method has a better transfer function than that of the constant step size 
RDC method. 
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Fig. 7.23. Diagram of the variable step size RDC system. 
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Fig. 7.24. The LPF causal impulse response. 
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when the variable step size is employed. 
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Fig. 7.26 The staircase representation of the causal impulse response of the LPF 

when the constant step size is employed. 
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Fig. 7.27. Transfer functions of constant and variable step size RDp LPF’s. 
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Fig. 7.28 A band suppression filter transfer function H n (f ) where 

BnU) = W) - H b (f ) 
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Fig. 7 . 29 . The impulse response of the band suppression filtejr . 
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Fig. 7.30. The staircase representation of the impulse response of the hand 
suppression filter when variable step size is employed . 1 
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Fig. 7.31. The staircase representation of the impulse response of the band 
suppression filter when constant step size is employed. 
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Fig. 7.32. Transfer functions of constant and variable 
step size RDC notch filters. 



VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The RDC system is a new realization of digital filters. The system can be 
used as a convolver or as a correlator. The RDC system has a simple layout for IC 
realization. When compared with other digital filter realizations, the RDC system 
appears to require less hardware and fewer connections. The system requires no 
synchronization, has no multipliers and is easily programmable. The RDC system 
provides real time operation (no off-line processing). There is no limitation on the 
time duration of the signal to be represented in matched filter applications. Simu- 
lation results indicate the RDC system achieves noise performance which surpasses 
that predicted by the theory of linear filters. 

B. RECOMMENDATIONS 

It is recommended that a RDC system be fabricated and tested to confirm 
simulation results and to provide additional insight into this method of achieving 
convolution. 

It is also recommended that additional study be done and the results be ex- 
tended to include adaptive delta modulation techniques. 

Also, studies of other noise reduction methods including additional work on 
slope filtering are recommended. 
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APPENDIX A 



IMPLEMENTATION OF A TRANSVERSAL FILTER 
USING DELTA CONVOLUTION 

Assume that we want to convolve the analog pulse signals h(t) and x{t) where 

h{t ) and x(t ) are the delta modulation versions of h(t) and x(t). Assume that h(t) 

and x(i) have number of bits n\ and n 2 respectively where 



n x < n 2 



Assume also that we have two registers Rh and R x where Rh holds the delta 
modulation bits of h{t), and R x holds the delta modulation bits of x{t) which have 
to be convolved with the contents of Rh ■ So, Rh has n x memory stages and R x has 
n 2 stages. 

By expanding Equation (3.18) we can write the result of convolution at time 
j -T a as 



y(t) ~ y(i • T s ) = c 1 



\ho 


-j-i j - 2 


j- 3 

+ E Xl + - 


j-n x 

..+ E xi 


{ 


.1=0 i= 0 


1=0 


i=0 




j- 2 


J— 3 


j —n i 


+ /ll 


E xi 


+ E i '+- 


• • + E x * 




i=0 


1=0 


1=0 






j- 3 


j-n l 


+ ^2 




E x « + - 
1=0 


. . + x *‘ 
i=0 



+ 



+ 

+ ^ni-l 



E X i 
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which can be rewritten as 



y(t ) ~ y(j ■ T 3 ) = Cl • < ho 



Xj—i + 2xj—2 + 3xj—o + 



j-n i-l 

+ ^1 • Xj-ni n i * y ^ 

2=0 



+h\ 



Xj- 2 + 2Xj_3 + 



+ (r»! - 1) • Xj- ni + (n x - 1) • ^ x i 

i=0 



+ 



+ ■ 



+ ^n t -1 



j-ni-1 

J2 xi 

i= 0 

(A.1) 



The products needed to implement Equation (A.l) can be represented as 
shown in Fig. A.l. A schematic diagram of a system which provides y(j • T s ) 
according to Equation (A.l) is shown in Fig. A. 2 where f s is the sampling fre- 
quency. 

In Fig. A. 2 we assume that Rh register stages can have any of three states 
— 1,0, or +1 and that it is bidirectional. Assume also that we have a logic which 
can do the following: 

a) During the time periods between even and odd sampling pulses the contents 
of register Rh move to the right. 

b) During the time periods between odd and even sampling pulses the contents 
of register Rh move to the left. 

Register Rh has 2ni stages. The first n\ stages have to be initialized with the 
delta modulation version of h(t ) which can be either +1 or —1. The last ni stages 
of register Rh have to be initialized with zeroes. 
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Fig. A.l. Diagram showing the products needed for delta convolution. 
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Fig. A. 2. Block diagram of the direct implementation of 
delta convolution. 



From previous assumptions, during each time period between even and odd 
sampling pulses the contents of Rh move to the right n\ times producing all prod- 
ucts needed to find the convolution, and during each time period between odd 
and even sampling pulses the contents of Rh move to the left rzi times producing 
all products needed to find the convolution. The contents of the accumulator are 
multiplied by the constant C\ to produce the convolution. 

Some of the disadvantages of this direct implementation of the delta convolu- 
tion are 

1 ) By looking at the columns of the XNOR gates of Fig. A. 2 , we see that 
the needed number n 3 of binary multipliers (XNOR) gates) is 



n 3 — n 2 + (^2 — 1) + (^2 — 2) + 

n 2 ri2 — Tii 

-5> X ■ 



• + (jl2 — Tli + 1 ) 



t=l 



i=l 



„ . "i K ) 2 

-m-nj + y- — 



2) The memory required is 2ni + ri 2 bits. 

3) -The high clock rate n\ • f s needed to clock Rh. 

It is clear that the required number of XNOR gates and the required memory 
will be large, especially when is large (which is a practical assumption). 
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APPENDIX B 



MODIFIED AND REDUCED DELTA CONVOLUTION METHODS 
1. Modified Delta Convolution 

We will derive the change in the discrete convolution A conv(i ) when the input 
register R x is fully occupied where 



Aconv(i) A conv(i ■ T a ) — conv I [i — 1] • T s 



(B.l) 



From Equation (3.14), at time t 



y(t) « y(k ■ T a ) = y[ [j + n x ] • T a ) = C x 



ho • [xj+ ni _i + 2xj+ Ul _2 4- 3xj+ ni _3 + 4- (ni — 2) • Xj + 2 4- (ni — 1) 

j 

Xj + 1 "I” Uj • ^ ^ 



+ h 2 • [ 



+ 



i=0 



+ h\ • £ xj+ ni _2 “I - 2 xj+ ni _3 4~ • • • 4- (nj 3) • Xj +2 4" (*^1 2) 



x j + 1 4- (nj - 1) ^2 x i\ 



i=0 



xj+m - 3 4- ... 4- («! - 4) • x j+2 + (ni - 3) 



xj -\- 1 4" (nx 2) i] 



i=0 



+ 

4- h ni _i • [ 
where 



£**•] 



!=0 



k = integer value of 



t 

Z 



= j+ ni 
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and where the needed products can be demonstrated as given in Fig. B.l. Collect- 
ing terms gives 



y(k • T a ) — y ([?' + ni] • T s ) — C\ 



Or 



ho • 

4- ^1 • 

+ ^2 • 
+ 



j +n i — 1 

53 Xi-(ni-i+j) 

i=j 



J-l 

+ ho - Til • 7> 
i= 0 



j + n x -2 

53 Xi ■ (ni - i + j - 1) 
i=j 

j+ni-3 

53 xi • Oi - * + j - 2) 

i=j 



j~ 1 

+ hi -(«1 - !) • 53 z * 
1=0 

i-i 

+ fr 2 (rci -2) • 53** 

«=o 



+ 



+ ^n!-l * 



j-l 



y; x » ■ ( n i - * + i - ( n i - 1)) + h ni -i • xi 






i= 0 



-t/(fc • T,) = C, • | £ (h 

' m=0 \ 

+ (X>) ■ (E **•[».- £! 

V.= 0 / V L=0 

At time (& + 1) • T s , 



j-fn x — m— 1 

53 Xi ■ (ni - i + j - m) 



( 5 . 2 ) 
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Fig. B.l. Diagram showing the needed products to find the convolution 

at time (j + n\)T 3 . 
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y ([^ + 1] ' T a ) — y ([ni + j + 1] • T s ) — C\ • ^ 

ho [•£ j+ni ”1" 2x j+m —l "H j^. ni _2 "I - H - (^i ~) ' xj + 3 

j + 1 ‘ 

+(ni — 1) • Xj+2 + • 53 Xi 



i=0 



+ hi • [ 



*^j+ n i— 1 ”h 2x j_|_ ni — 2 + • • • + ( n l — 3) • Xj+ 3 

7+1 ■ 

+(ni — 2) • Xj+2 + ( n i ~ 1) • 53 x i 



x=0 



+ ^2 • [ 



Xj-f m —2 + ... + (r»i - 4) • Xj +3 

j+i 

+(ni - 3) • x J+2 + (ni - 2) • 53 x t - 



i=0 



+ 

+ h ni _i • 

Therefore;, 

y([fc + 1] -T,) = y ([nj + j + 1] • T s ) = C\ 
ho • 



i+i 

^Xi 

t=0 



7+ n i 

^2 x * ■ ( n i -i+j + 1) 
«=j 



+ h 



i * 



7 + n l -1 

53 Xi-(ni-i+j) 

*=i 



+ ho I — x j +ni ■ y> 

t=0 

7-1 

+ hi I Xj + [n x - 1] • 53 Xi 



i=0 



+ 
+ 



7+ n i “2 

53 H • (n! - i + j - 1) 

i=j 



7-1 



+ h 2 I -Xj + [ni - 2] • 53 x,' 



1=0 



/ 


*7 + 1 


1 hni —l 


Xi • (n i - i + j - [ni - 2]) 

_* = j 



+ h n . _i 



i-i 

— xy + 53 x,- 
i=0 J 



(5.3) 



Or, 
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v([t + l]-T.)-C,-( J 2 l A » 



L m=0 



j+n i — m 

x; • (n 1 - i + j - m + 1) 
i=i 



O'-i 



.i=0 



"ru-l 

L=0 



Til — 1 



+ X Xi ' X h L * hi - £] ) - Xj • Yj h i f • ( BA ) 



i=0 



From Equations (B.2) and (B.4), now form A conv(j + nj +1) where - 



A conv(j + ni + 1) = conv ([j 4- n\ + 1] • T s ) — conv (\j + n\] • T s ) 

j + n i — to 

X x » ■ ( n i — * + j — m + 1) 

*=> 



^ * ( { X ( 

m=0 



n t -l 

— Xj • A; 

i=0 



{m-i 

X i hm 

m= 0 



7 + ni — m —1 

X • (wi - i+j-m ) 
i=j 




= Ci 



Til — 1 

X 

m=0 



}-+!>) 



= C 2 



j+n x — m 

X Zi • (ni - i + j — m + 1) 

»=? 

j+711— m — 1 

- X * («i - * +j - m) 

»=j 

n t — 1 f / i+ni — m— 1 

X 1 ( X x »" K ni - * + j - m + !) - (Til - * + j 

< \ i=j 
m-i 

— x j ■ X 

i=0 

ni- 1 / j + Ti! — m — 1 \ ni- 1 

^ ] A m | ^ ^ 3 Tj + I — Xj • ^ ^ 

i=j / i=0 



_ m=0 



c x - 

"I” X j+n x — 



-m)] 



m=0 



= C X 



ni — 1 /i+ni — m 

XM X xi 

m = 0 \ isj 



ni— 1 

” h* 

i=0 



ni- 1 /j+ni— m 

^ * ( X ( X x » 

m=0 V 1 = 7+1 



(B.5) 
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Or, 



A conv (j + ni + 1) = Ci • 




< Rh(m ) > • 



ni —1 







(5.6) 



where < Rk(i) > means the contents of cell i of register k where k can be h or x. 

We can see that the change in the discrete convolution is dependent only on 
the last n-y samples of the input and is independent of all other samples. Also we 
can see that all the weights equal unity. The needed binary products to find the 
change in the discrete convolution are represented in Fig. B.2. 



2. Implementation of a Transversal Filter Using the Modified Delta 

Convolution Method 

To implement a transversal filter using Modified Delta Convolution, two 
methods can be used. 

a. The First Method 

In the first method (Fig. B.3), the needed binary multipliers (XNOR 

gates) = ny + (n{ — 1) + («i — 2) + + 1= n d 2 +( n i - 1 ) = m (in + 1 ) anc J the 

needed binary storage =2 ny . 

b. The Second Method 

In the second method (Fig. B.4) we will use the fact that we need to 
perform ny binary multiplications only and to store all the needed binary multi- 
plications which were performed in the preceding steps. In this case the needed 
number of binary multipliers = n y , and the needed storage is as follows: 

2ny + (ny — 1) + (ni — 2) + + 1 = n t + n 1 ■ "-, 1 + - D = "i L tti±d l . Despite the fact 

that in the second method the needed hardware is reduced, we still have to use a 
number of XNOR gates or a number of binary storage locations proportional to 
(ni) 2 . Because of this we introduce the Reduced Delta Convolution method. 
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Fig. B.2. Diagram showing the needed products to find the change 

in the discrete convolution. 
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Fig. B.3. Diagram showing the first method of implementing 
modified delta convolution. 
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Fig. B.4. Diagram showing the second method of implementing 
modified delta convolution. 



3. The Reduced Delta Convolution (RDC) system 

In this section we introduce the reduced delta convolution (RDC) method 
where the needed number of binary multipliers equals the filter length n x and the 
needed binary storage equals 2 n x 4- 1. 

Let us define A2 conv(k ■ T a ) as follows 

A2 conv(k • T 3 ) = Aconv(k • T 3 ) — A conv ([Ar — 1] • T 3 ) 



where A conv(k ■ T a ) is the change in the discrete convolution which is defined in 
Equation (B.l). 

We now find A2 conv(k • T 3 ) when the input register is fully occupied. At 
time t where 



k = integer value of 




= j + n i 



we know that 



A conv(k • T 3 ) = A conv([j + ni] • T a ) = C\ 
ho [^j+nt-l + Xj+nx- 2 + Xj+ ni - 3 + + Xj- 1 + X j] 



h\ [ _2 “I" 3 4" 4" Xj — 1 4“ Xj] 

+ h 2 [ Xj +Tll _3 4- 4- Xj - 1 + Xj] 



+ 

+ 

4- Aim-i [ 




(B.S) 
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and 



A conv([k + 1] • T a ) = A conv{\j -f ri\ + 1] • T a ) = C\ • 

ho [xj. j. ni d- ij+m-i d - ^7+ni —2 d - d - *7' d- *7+1] 

d-^l [ *7+ni-l d- Xj+ ni -2 d- d- d- £7 + 1] 

d"^2 [ x j+ni ~2 d* d- £7 d* £7+1] 

+ 



d- 

d-^m -1 [ 



*7+1 ] 



(B.9) 



From Equations (B.7), (B.8), and (B.9) we have 

A2 conv ([& + 1] • T s ) = A conv ([& + 1] • T s ) — Aconv ( k • T s ) 

A2 conv (| j + ni + 1] • T s ) = Aconv ([7 + n\ + 1] • T s ) — Aconv ([? d- «i] • T 3 ) 

— C*! ' \ -ho Xj- j-fij d* h\ Xj +Tll _ j h 2 Xj^-rti —2 d“ • • • 



d- h ni _i Xj+i Xj 



n l — 1 

E h < 

. i=0 



(B.10) 



which can be represented by Fig. B.5. Referring to Fig. B.5 and Equation (B.10), 
then 

( r»i — 1 

A2 conv([j + ni + 1] • T a ) =C X • < ^ < R h ( i) > ■ < R x {i ) > 



t i=o 



n \ -1 



< R x (p, 1) > • ^ ^ h{ 



(B. 11) 



i=0 



From Equation (B.l), 



Aconv ([j + ni + 1] • T s ) = conv([j + nj] • T s ) + A conv(\j + n x + 1] • T s ) (£.12) 
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and from Equation (B.7), 



Aconv([j -f ni + l] • T a ) = A conv([j + «i] • T a ) + A2 conv([j + m + 1] • T a ) 

(B. 13) 

where A'2conv([j + n\ + 1] • T,) is given by Equation (B.ll) and where at 
time t = 0, conv{ 0) — 0. 

Equations (B.ll), (B.12), and (B.13) define what we call the reduced delta 
convolution (RDC) system. See Fig. 4.2. 
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Fig. B.5. Diagram showing the creation of A2 conv (j [j + nj + 1] • T 3 ). 
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APPENDIX C 



A CONVENTIONAL MATCHED FILTER 



Fig. C.l shows a block diagram of a matched filter where 



y(t)= / h(r ) • w(t — r)dr 



-J 

-J 



x(Td — t ) • w(t — r)dr 



(C.l) 



where x(t ) is the signal to which the matched filter is matched and Tj is its time 
duration. 

The first approximation involves writing a sampled-data form of Equation 
(C.l). See Equation (3.15). 



L- 1 

Vk ~ T s • 'y ^ i * W k—i 

t— 0 



(C.2) . 



where 



y ; - = y(j • T s ), wj = w(j • T 3 ), and xj 4 x(j ■ T 3 ) 



The sampling rate = jr must exceed the Nyquist rate. L is the number of samples 
and so Tj in Equation (C.l) equals L ■ T s . 

The samples xj and Wj axe now represented by binary words having elements 
1 or 0 identified by x ™ and wj where 

OO OO 

x, = E x T ■ 2 "”- = E “ l ■ 2 '" ( c - 3 ) 

m=0 n=0 

where x J 1 and w J G {0, 1 } . 
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w( t ) 






h( r ) = 

X( T d - T) 



y( t ) 



Fig. C.l. Block diagram of a matched filter. 
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We now approximate the sum in Equation (C.2) by truncating xj to M digits 
and Wj to N digits. As a result of this 



A /-1 



N - 1 



L— 1 



« « t, • y , E 2 ~ (m+n) • E ( c - 4 ) 

m= 0 Ln=0 i=0 

Correlating the m— digits of the binary words of the samples of x(-) with 
the n— digits of the binary words of the samples of w(-) gives the digital corre- 
lator (m, n) given in Fig. C.2. At each sample time the most recent sample’s 
n— digit enters the left most stage of the n— on-line shift register. The m— dig- 
its of the sample xj are permanently stored in the stages of the reference shift 
register. But, each sample Wj has M bits and each sample xj has N bits. So, 
to obtain yjt we have to perform M ■ N digital correlations, where the results 
of correlation are weighted by the factor T 3 ■ 2“( m+n ). We use one reference 
shift register, which has M • L binary stages, for all the M ■ N digital correla- 
tors. 



Fig. C.3 shows a block diagram of the DMF for arbitrary N and L but with 
M equals to 1 because a block diagram for arbitrary M, N and L showing all inter- 
connections is cumbersome. In Fig. C.3, for convenience, the gain factor T a called 
for in Equation (5.8) has been left out of the weights 2~( m+n \ so the output in 
fact is . 
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the nth onjine shift register 







Fig. C.2. Digital correlator (m, n). 
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Fig. C.3. Block diagram of a DMF for arbitrary N and L, where M — 1. 
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APPENDIX D 



DC OFFSET LEVEL 

After the input expires, the zero voltage input to the delta modulator will be 
represented by a successive positive and negative pulses about zero. From Equation 
(3.18) the output of the convolver equals 



ni-l ( 

y(t ) « y(j ■ T s ) = Ax ■ Ah ■ T s ■ ^ < 



k=o K Lp=o 






j-k-i 



X i 



L i=0 



! j — 1, 2, 3, . . . 



But the input signal delta modulation pulses X{ have successive positive and neg- 
ative values. Therefore, 



j-k-i 



— A. — 1 s 

E*-{ 

i=0 



or — 1 for j — k = odd 
for j — k = even 



(D.l) 



From Equation 3.11 



h(T)Kh 3 (r) = Ah-Y / h p; k ■ T s < r < (k + 1) ■ T 3 (D.2) 

p = o 

Therefore, the output level at a given sample time j ■ T 3 when j is odd (k even) 
becomes 



\ 

E 

fc=° / 

\fc even / 



(±A;r) 



At T 3 seconds later, j will be even ( k odd) and so the dc offset value is 



ni -1 



T 3 - £ h(k-T s ) .(±Ax) 



k = 0 

^ k odd 
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Because T a is small when we use DM , then we can expect that the distortion 
will be constant and approximately equal to 



±- 



Ax 



(total area under h(t)) 



(D. 3 ) 



„ n , -1 

since T 3 • ^ h(k ■ T 3 ) is approximately half the area of h{t). 

k = 0 

k odd or even 
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